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ABSTRACT 

The  purpose  of  the  present  work  is  to  study  surface 
waves  in  rectangular  channels  of  finite  depth.   In  all  cases 
we  consider  two  dimensional  flow  of  an  incompressible  and 
nonviscous  fluid.   We  also  suppose  that  there  is  surface 
tension  present  in  it. 

In  the  first  chapter  we  discuss  the  feasibility  of  a 
steady  flow  consistent  with  the  linearized  problem,  and  study 
the  corresponding  radiation  conditions.   Such  a  flow  is  found 
to  be  possible  in  both  the  supercritical  and  the  subcritical 
cases  but  not  in  the  critical  case.   In  the  supercritical  case 
the  disturbance  dies  both  upstream  and  downstream;  in  the 
subcritical  case,  in  general,  it  does  not  vanish  downstream. 
The  second  chapter  deals  with  a  study  of  solitary  and 
periodic  waves  which  result  when  a  stream  of  uniform  flow  is 
disturbed  by  an  arbitrary  variable  pressure  on  its  free  surface 
We  obtain  a  general  differential  equation  of  the  motion. 
Explicit  wave  solutions  are  obtained  by  considering  two  special 
pressure  distributions. 

In  the  final  chapter  we  study  the  progressive  wave 
solutions  of  the  zero  order  shallow  water  equations.   In  the 
subcritical  case,  we  obtain  a  wave  with  a  single  hump  and  a 
single  hollow  on  the  two  sides  of  the  point  of  disturbance. 
In  the  critical  and  the  supercritical  cases,  we  obtain  periodic 
waves. 


Ill 


Introduction 

The  purpose  of  this  paper  is  to  study  surface  waves  in 
an  infinite  channel  of  finite  depth.   The  channel  in  all  cases 
is  assumed  to  have  a  horizontal  bottom  and  parallel  walls. 

The  phenomenon  of  surface  waves  is  approached  through 
two  approximations,  i.e.  the  amplitude  of  the  wave  is  small 
which  yields  a  linear  theory  or  that  the  depth  of  the  water  is 
small  compared  to  some  other  significant  length  such  as  the 
radius  of  curvature  of  the  water  surface  which  results  into 
a  non-linear  theory.  In  what  follows,  the  first  chapter  is 
based  on  the  former  and  the  remaining  two  chapters  on  the 
latter. 

In  the  first  chapter  we  consider  the  feasibility 
of   the         propagation  of  disturbance  as  steady  waves. 
The  discussion  is  carried  out  within  the  framework  of  linear 
theory.   It  is  assumed  that  surface  tension  is  present  in  the 
water.   In  cases  where  such  waves  exist  we  seek  the  corre- 
sponding radiation  conditions. 

This  question  becomes  relevant  because  as  Stoker  [ 0  ] 
points  out,  the  steady  state  formulations  are  not  natural  in 
Newtonian  mechanics.   The  unnaturalness  stems  from  the  fact 
that  whereas  Newtonian  mechanics  deals  with  the  initial  value 
problems,  the  steady  state  problems  are  those  in  which  time 
does  not  occur  at  all.   The  assumption  underlying  the  latter 
formulations  is  that  after  a  long  period  of  time,  the 


propagation  of  disturbance  tends  to  be  steady.   The  position 
taken  In  this  chapter  Is  that  It  Is  not  a  priori  obvious  and 
therefore  It  should  be  Investigated.   Indeed,  our  results 
reveal  that  when  IT  =  gh,  the  steady  state  does  not  exist. 
It  exists  only  for  the  case  U  >  gh  and  U  <  gh. 

This  discussion  gives  rise  to  the  question  as  to  how 
far  the  steady  state  problems  produce  a  reasonable  solution. 
It  Is  a  matter  of  common  experience  that  such  formulations 
In  an  Infinite  domain  require  more  stringent  conditions  than 
those  In  the  corresponding  time  dependent  problems.   These 
are  the  so-called  radiation  conditions.   Aside  from  an 
Intelligent  guess,  there  Is  no  obvious  way  to  find  these  ex- 
cept by  solving  the  problem  Itself.   In  the  present  Investiga- 
tions these  have  been  found  by  actually  solving  the  Initial 
value  problem  and  passing  to  the  limit  as  t  — >  qd  .   We  find 
that  the  disturbance  must  die  out  upstream  and  downstream  In 
case  IT  >   gh.   When  IT  <   gh.  It  vanishes  upstream  but.  In 
general,  not  downstream. 

These  Ideas  were  originated  by  Stoker  and  he  solved  the 
above  problem  [9].   The  present  Is  a  generalization  of  the 
same  when  water  admits  surface  tension. 

The  second  approximation  as  alluded  to  above  develops 
the  classical  long  wave  theory.   As  propounded  by  Lamb  [  4  ] 
It  has  Its  basis  In  hydrostatic  pressure  law.   In  Its  modern 
version.  It  was  obtained  by  Frledrlchs  In  an  appendix  to  a 
paper  by  Stoker  [  2  ],  as  lowest  order  approximation  In  the 


shallow  wave  theory.  We  follow  a  similar  procedure  to  dis- 
cuss solitary  and  cnoldal  waves.  This  discussion  forms  the 
contents  of  the  second  chapter. 

A  solitary  wave  is  a  wave  having  a  symmetric  form  with 
a  single  hump  and  which  propagates  at  uniform  velocity  with- 
out change  of  form.   Cnoldal  waves  are  periodic  waves.   The 
solitary  wave  was  first  observed  by  Russell  in  l844  [5].   The 
mathematical  discussion  of  it,  i.e.  its  form  and  the  velocity 
at  which  it  propagates,  was  initiated  by  Boussinesq  [1].   He 
was  followed  by  Rayleigh  and  several  other  authors.   In  its 
next  stage  of  development,  Kortweg  and  DeVries  discussed  the 
general  class  of  waves  which  propagate  without  change  of  form. 
Their  investigations  produced  both  the  solitary  wave  and  the 
cnoldal  wave. 

The  basis  of  the  development  here  is  the  shallow  water 
theory  as  given  in  Stoker  [2].   (See  also  Stoker  [6].)   The 
method  essentially  consists  in  making  a  stretching  trans- 
formation of  the  variables  occurring  in  the  problem,  assuming 
asymptotic  expansions  of  the  unknowns  in  terms  of  a  small 
parameter  and  solving  for  the  functions  appearing  in  these 
expansions.   Keller  used  this  technique  to  discuss  solitary 
and  cnoldal  waves  [?]•   A  variant  of  it  was  used  by  Stoker 
and  Peters  to  discuss  solitary  waves  in  a  medium  of  varying 
density   [8].     Since  then  it  has  been  studied  in  numerous 
papers.   A  simplification  Introduced  in  all  these  cases  is 
that  the  pressure  on  the  free  surface  is  taken  to  be  zero. 


Instead,  we  consider  here  that  an  arbitrary  pressure  Is 
applied  at  the  surface.   The  result  Is  a  differential 
equation  of  Kortweg  and  DeVrles  type  with  a  variable 
pressure  term  In  It.   Explicit  solutions  of  this  equation 
are  obtained  In  terms  of  solitary  and  periodic  waves  by 
considering  two  particular  pressure  distributions. 

The  shallow  water  equations  In  the  lowest  order  of 
approximation  are  known  to  have  the  uniform  flow  as  the 
only  smooth  solution.   In  the  final  chapter,  we  discuss 
the  nontrivlal  progressive  wave  solutions  when  surface 
tension  is  present  in  water.   The  discussion  involves 
a  parameter  (uQ-a)/Y/gh  depending  on  the  original 
unperturbed  flow  and  the  depth  of  the  channel.  When 
this  parameter  is  less  than  one,  we  get  a  wave  starting 
from  the  point  of  disturbance  and  spreading  over  the 
whole  right  half  and  left  half  axis  with  respectively 
a  single  hump  and  a  single  hollow.   In  the  critical  as 
well  as  in  the  supercritical  case,  i.e.  (uQ-a)  /gh  _>  1, 
the  wave  becomes  a  periodic  one. 


Chapter  I 
Linear  Theory  of  Steady  Waves  In  Channels  of  Finite  Depth 

Formulation  of  the  problem 

We  consider  two  dimensional  flow  of  an  Invlscld  In- 
compressible liquid  In  an  Infinitely  long  channel  of  finite 
depth,  h.   The  liquid  has  surface  tension  T.   We  suppose  that 
It  Is  Initially  moving  with  velocity  U  parallel  to  the  sides 
of  the  channel  until  some  disturbance  Is  created  at  the  free 
surface.   The  problem  Is  to  study  the  feasibility 
of  the     propagation  of  this  disturbance  as  a  steady  wave. 

In  order  to  obtain  a  reasonable  solution  of  a  problem 
posed  In  an  Infinite  range.  It  Is  necessary  to  Impose  radia- 
tion conditions  at  ±cx)  .   Normally,  In  a  physical  formulation, 
boundedness  conditions  of  some  kind  should  suffice.   It  Is 
common  experience,  on  the  other  hand,  that  In  time-Independent 
cases,  more  stringent  conditions  are  required.   As  to  what 
these  conditions  should  be  Is  not  a  priori  obvious.   We  seek 
to  find  these  radiation  conditions  for  the  problem  mentioned 
above  In  the  event  the  steady  state  does  exist. 

Such  an  Investigation  has  already  been  carried  out  by 
Stoker  [  9  ]  for  the  case  when  the  fluid  has  the  surface 
tension  0.   The  present  Investigation  is  a  generalization  to 
the  case  of  nonzero  surface  tension.   We  study  the  linearized 
problem. 


We  choose  our  x-axls  parallel  to  the  sides  of  the 
channel  In  the  equilibrium  plane,  the  y-axis  is  chosen 
vertically  upward.   We  assume  that  the  disturbance  is 
created  by  a  pressure  7T(x,t)  on  the  free  surface. 
Since  the  motion  is  originally  irrotational,  it  will 
remain  so  always.   Hence  there  exists  a  velocity 
potential  of  the  motion.   We  denote  by  J,(l)  the  complete 
and  the  disturbed  flow  respectively.   Thus  J  =  ^   +  Ux. 
Let  r[{x,t)    be  the  free  surface.   The  equations  of 
motion  are: 

A  ?(x,y;t)  -  0 


Pq  1^2 

—  +  $1-  +  gT  +  g^ix  ^^y^  "  const,   at  y  =  rj 


where  p^  is  the  pressure  at  the  surface  of  the  liquid. 

Tri 
'xx 


'X 

Eliminating  p^^  between  the  two  foregoing  equations  we  have 


The  kinematical  condition  at  the  surface  is 


^t  ^  ^x^x  "  ^y  "  °  '      at  y  =  ri  ; 
^  =  0       ,     at  y  =  -h. 


The  assumption  of  linearization  physically  implies 
that  the  disturbance  generated  is  so  small  that  the  terms 
of  order  higher  than  the  first  are  neglected.   The 
linearized  equations  are 

(1)  A(i)  =  0 


(2)  p  +  S"!  "^  'i't  "^  ^'^'x  "^  I  "^xx  =  °      at  y  =  0 


where  the  constant  on  the  right  is  absorbed  in  ^ 


(3)  T)^  +  U^x  =  "^y  at  y  =  0 


(4)  "^y  =  °  at  y  =  -h, 


At  time  t  =  0,  we  assume 

<[)(x,y,0)  =  Ti(x,0)  =TT(x,0)  =  0 

whence  from  (2),    assuming  that  it  holds  on  the  free  surface 

([)^(x,y;0)  =  0 
Prom  now  on  we  assume 

TT(x,t)  =  TT(x)  for  t  >  0. 

Finally,  the  conditions  at  +oo  for  ^   are  considered  in 
the  existence  of  the  Fourier  Transform  (j)  of  (J)  and  its 
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derivatives.   The  conditions  on  TT(x,t)  are  also 
prescribed  through  its  Fourier  transform  and  are  kept 
open  for  the  time  being. 


®  -Isx 


(j)  =  — ^   /  (l)(x,y;t)  e    dx 
/27r 


-co 

00 


■isx 


~ /  ^^(x,y,t)  e    dx 


1/2?  s  _^   ^ 

00 

^  =  -^   /"(})   (x,y,t)  e    dx. 


-ISX 

"  '    '^  -^  '"" 


Also 


®  -isx 


(j)   =  — ^   /  h(x,Y,t)    e  dx, 

/     -oo 

Thus  (1)  becomes 


(5)  ^yy  -  ^  '^^  =  °  • 


Sjjnllarly    (2),    (5)    become 

(6)  —-  +  g^  +  ^t  +  u(i(t)s)   -  s^  -^  ^  =  0 


(7)  Tit  ^  ^^^^   "  *^y  "  ° 

where  t^",]!  are  the  transforms  of  Ti^fT  ^esps^^^v^^y  ■ 
Solving  equation  (s),  we  get 
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(8)  <t>Cs,y;t)    =  A(s,t)    cosh   s    (y+h) 

Ti(s^  ^  -g)    =—  +~i^   +  Is'Ju 


Substituting   for  T],r\      In   equation    (7)    we  get 


(9)        ^^^   +  21sU(t)     -  sVct,  -   ^^  fs^  I  -g)   =  -   isu  n   . 

y        P  p 


Substituting  for  ^   from  equation (8),  and  evaluating  at 
y  =  0,  we  have: 


A^^cosh  sh  +  2isUA^  cosh  sh  -  A{s^U^cosh  sh  + 

+   s(s  -  -g)slnh  sh}  =   -IsU  11  . 
P  p 

^tt   +  2^^U^t   -  A(sV+  s(s2  1  -g)tanh  sh]    =      -isU  Asech  sh 


In  order  to  solve  the  homogeneous  equation,  we  have 
to  solve  the  quadratic: 

D^  +  21sUD  -  (sV+sCs^l  -g)tanh  sh]  =  0  . 
D  =   -isu  +  {sfs^l  -g)  tanh  sh}^/^  _ 

r 

The  solution  of  the  equation,  therefore,  is 


A(s,t) 


isuTT 

p  cosh  sh 


2„2   ,    2  Tn, 
s  U  -s(g-s  — )tans  sh 


+  B  e 


-lt(sU+  /s(g-s^  -)tanh  sh) 

r 


-lt(sU-  /s(g-s^  -)tanh  sh) 
+  C  e  ,  for  |s|<  /■£§  ; 


'- V  T 


and 
A(s,t) 


IsUTT    

p  cosh  sh  ^  ^2„2.  /  2  T  „^^^^,     „>. 
r  s  U  +s ( s  —  -gjtanh  sh 


< 


(-lsU+  /s(s^  -  -g)tanh  sh)t 


+  B,  e 


+  C^  e 


(-isU-  ys(s^  -  -g)tanh  sh)t 


,for  |s|>7-^ 


B,C,B-,C,  are  constants  of  Integration  to  he  evaluated 
from  the  initial  conditions  A(s,0)  =  0  and  A^(s,0)  =  0 
These  follow  from  equation  (8)  and  the  fact  that 
(|)(s,0)  =  (i)^(s,0)  -  0.   Thus 


B  +  C  =  - 


p  cosh  sh 


i 


2  2       2  T 
s  U  -s(g-s  — )tanh  sh 


B(sU  +  l/s(g-s  -)tanh  sh  )  +  C(sU-  ^s(g-s  ^)tanh  sh  )  =  0 


c  = 


-isuTT 


1 


/     2  T  I     /  2  T 

2p  cosh  shysCg-s  — )tanh  sh  sU-  /s(g-s  — )tanh  sh 

P  ^  P 
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B  - 


isuTT 


2p    cosh  sh 


ys(g-s^  -)tanh  sh 


sU+  /s(g-s^  -)tanh  sh 


Similarly, 


B. 


-  isuTT 


2p    cosh  shys(s     -  -g)    tanh  sh 

r 


r 


IsU  +  /s(s^  -^  -g)tanh  sh 


^1  = 


2p  cosh  sh/s(s  -^  -g)tanh  sh 

r 


^   -isU  +  /s(s^  -  -g)tanh  sh 

2  2       2  T 
^  s  U   +  s(s  -^  -g)tanh  sh 


> 
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isuTT 

p  cosh  sh 


s  U  -  s(g-s  — )tanh  sh 


-i{sU+  ys(g-s^  -)tanh  sh  }t 
+  ^ 


2  T  Z'  2  T 

2  /(g-s  — )s  tanh  sh  { sU  +  /s(g-s  — )tanh  sh  } 

r  F 

-l{sU-  >4(g-s^  -)tanh  sh  ]t 


2  /(g-s^  -)s  tanh  sh  { sU-  /s(g-s^  -)tanh  sh} 

r  r 

for  I  s  I  ^  7-^  =  a 


A  =   and 


p  cosh  sh 


s  U  +s(s  —  -g)tanh  sh 

r 


(lsU+  /sfs'"^  -  -g)tanh  sh)  t 


2  /s(s^  -  -g)tanh  sh  {-isU+  /s(s^  -^  -g)  tanh  sh  } 

r  r 

(-IsU-  /s(s^  -  -g)tanh  sh  )t 
e__ ^ 

2  /s(s^  -  -g)  tanh  sh  [  lsU+  /s(s^  -  -g)tanh  sh  ]_ 

r  r 

for  I  s  I  >  y-^  . 
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We  can  solve  for  (j)  by  Inverse  transformation 

00 

1       P  ISX 

<l)(x,y;t)  =  /  A(s,t)  cosh  s(y+h)  e    ds 

The  convergence  of  the  above  Integral  for  a  finite  t 


<t,(x,y,t)  =-^   r   l3UTTco3h  s  (y+h)  e^^^ 
^/^    J  P  cosh  sh 


/2i    ^_ 


a 


s%^-s(g-s^  ■5^)tanh  sh        2  /s(g-s^  -)tanh  sh 


/' 


-l{sU+Xs(g-s^-)fcanh   sh}t  -1{  sU-/s  (g-s^ --)  tanh  sh}  t 

P  c.  P 


,  -mzz^mzzzzzz^i — r 

2   T  /  p    fp 

sU   +ys(g-s     — )tanh  sh  sU-  /s(g-s     — )tanh   sh       J 


/  -a      oov         _ 
I      1   {    r+    MlsUTTcosh   s    (y+h)    e^^^ 

Vp^^V  -'       ^    I  P   cosh  sh 

'       Voo      a/ 


ds 


+ 


_s^U^+s(s^  -^  -g)tanh  sh  2  /s(s^  -  -g)tanh  sh 


15 


-isUt+t/s(s^  -  -g)tanh  sh    -IsUt-t/s  (s^  -  -g)tanh  sh 
-  e ^ e ^ I 

-lsU+)4(s^  -  -g)tanh  sh      lsU+/s(s^  -^  -g)  tanh  sh 


=  (})  +  (t)    +  i  where  a  =  ^^  , 

^a   ^-00    ^oo  T 


ds 


where  ^     is  the  first  integral  from  -a  to  a  and 
a 

d)    and  i   are  the  remaining  two  integrals  in 

^-00        ^00 

(-00, -a)  and  (a,oo)  respectively. 

We  first  consider  (b  .   The  function  outside  the 

a 

parenthesis  does  not  have  any  singularity  within  the 
interval  of  integration.   We  claim  that  the  function 
inside  the  parenthesis  is  bounded.   This  can  be 
demonstrated  by  expanding  the  exponentials.   Thus 


+ 


j^^-s(g-s^  -)tanh  sh   2  /s(g-s^  ■^)tanh  sh 
P  P 


e 


-i(sU+/s(g-s^  -)tanh  sh)t   -it (sU-/s  (g-s^  |)tanh  sh 
P  p  ^ 


> 


sU  +  /s(g-s^  -)tanh  sh      sU  -  /s(g-s  -)tanh  sh 
°    p  p 


+ 


s^U^-s(g-s2  ^)tanh  sh   2  /(g-s^  -)tanh  sh 

P  r 
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(sU-  /s(g-s^  -^)tanh  sh)  { 1-lt  (sU+  /s  (g-s'^  -^Hanh  sh) 


2  t2(sU+/s(g-s2^)tanh  sh)^+.  .  .]-is\J+   /s  (g-s^  ^')  tanh  sh) 


'{1  -  it(sU-/s(g-s2|)tanh  sh)  -  1  t^CsU-  /s(g-s^  ^)tanh  sh)^ 


+  .  .  .} 


—1   2  2  ,         2  T 

s  U  -s(g-s  -)tanh  sh 


.2,^2 


2  T, 


1- 


s  U  -s(g-s  -)tanh  sh     ^  ^/~,        2~T^ 

P         L   2  /s(g-s''  -)tanh  sh 


•j2t  ys(g-s  |)tanh  sh  (l+t^Cs V-s(g-s^  ^)tanh  sh))  ' 


2..2 


2  T, 


=^  1  +  [s  U  -  s(g-s  ■i)tanh  sh}  +  ... 


That  is,  all  the  singularities  of  the  Integrand  are 
cancelled  out.   Since  s  takes  on  only  finite  values, 
the  function  in  the  parenthesis  is  bounded. 
Let  its  bound  be  M. 


M 


27r 


f     sU(Tr(s)+TT(-s))  cosh  s(y+h) 

J  p  cosh  sh 

0 


ds  , 


15 


which,  being  finite.  Implies  that  ^      converges. 

s. 

In  the  case  of  (|)   ,  the  singularity  at  the  point 
s  =  a  is  a  square  root  singularity  and  is  therefore 
integrable.   On  the  other  hand,  for  large  |s|. 


s  cosh  s  (y+h)  _  ^  e^(y+h)^^-s(y+h)       [^ 

i 1 —   O        _  1 \ ^^^   b   C 

cosh  sh         ^sh  ,   -sh 

e   +  e 


The  integrand  in  the  absolute  value  behaves  as 


sU 
)/27r 


ITT(s) 


.y|s 


s^U^+s(s^  -  -g)tanh  sh 


+ 


2  /s(s^  -  -g)tanh  sh 

r 


t  ys(s  -  -g)tanh  sh    -t  /s(s  -  -g)tanh  sh 
;        P  +  e         P 


/s^U^  +  s(s^  -  -g)  tahn  sh 


-TTr      \  \  y  S     t  S 

I  I  (s  )  I  se-^  '  '  e 


5/2 


for  large  s . 


The  integral  ^  then  will  be  convergent,  provided  TT(s) 


oo 


behaves  like  e    for  large  values  of  s .   Indeed  the 

subsequent  analysis  reveals  that  TT(s)  should  be  of  the 

2 

form  e  ^  g(s),  where  g(s)  remains  bounded  at  the  two 
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infinities.   Henceforth,  therefore,  we  assume  that 
_       _  2 
TT(s)  =  e   g(s),  where  g(s)  is  analytic  for  all  s 

and  is  bounded  at  +00.   The  function  g(s),  for  instance, 

can  be  a  periodic  function.   The  convergence  of  ^  can 

-OD 

be  discussed  the  same  way. 
We  write 


GO 


^^   =  -^Z-  f    A(s,t)s  cosh  s(y+h)e^^^  ds 


27r 


-OD 


This  Integral  then  gives  the  velocity  of  propagation 
of  the  disturbance  at  any  finite  time  t.   For  the 
existence  of  steady  waves  we  have  to  study  (J)   after 
a  long  time  over  the  whole  of  the  x-axis   i.e.lim  (J) 


lim  ^ 
t-^co  ^ 


lim 

27r  t-^00 


a 


t->aD 


isx 


A(s,t)  s  cosh  s(y+h)e    ds 


A(s,t)  s  cosh  s  (i'+h) 


xsx 


ds 


-/2^ 


\     +  ^-00  -*-  ^OD 


I„  - 


lim   ris^uTTe^^^  cosh  s(y+h) 

t-*OD  ^ 

-a 


P  cosh  sh 


1 


2„2    T,  2   2^,   ,   , 
s  U  -s  —fa  -s  Jtanh  sh 
P 


+ 


-if  t    -if_t 
e        e 


/T     2   2 
—  s(a  -s  )tanh  sh 

r 


f 


+ 


ds  , 
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where 


f_^=   sU  +  V^(g-s^  -)tanh  sh 

r 


+ 


(11) 

Let 

Independent  part 


•'■a  "  ■'■as  ^  -"-at  '  ^^^^^   ^as  ^^  ^^^   ^^^ 


We  first  consider  I  , . 

at 

This  requires  the  Investigation  of  the  zeros  of  f , j f 
We  write 


F(s)  =  f ^  •  f . 


s^U^  -  (g-s^  -)  s  tanh  sh 

VO         p 


=  sg 


2  2 

^  hs  -  (1-  -^-^)tanh  sh) 
gh  gp 


F(s)  =0,  If  s  =  0,  or 


U^   V,      /^n     2  T  \  ,   , 
— T-  •  hs   =  I  1  -  s  —  tanh 
gh         V        gp  / 


hs 


It  Is  obvious  that  this  equation  also  has  a  zero  at  s  =  0 , 
Hence  P(s)  has  double  zero  at  s  =  0 .   For  s  7^  0,  we  write 


U 
gh 


hs  coth  hs  = 


1  -  s  h  


y 


U 

gh 


pgh 
say. 
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Thus  the  zeros  of  P(s)  other  than  s  =  0  are  the 
intersection  of  the  curves: 

y  =  hs  coth  hs        (^) 


^  pgh 


Li^) 


It  follows  that  for 


>  1,   F(s)  has  a  zero  of  second  order  at  s  =  0 


gh 

u2 

—r  =   1,      F(s)  has  a  zero  of  fourth  order  at  s  =  0 
gh 

u2 

— T-  <  1,   F(s)  has  a  zero  of  second  order  at  s  =  0 

and  two  simple  zeros  at  s  =  ^,-^;     Ip  I  <   pg/T.   We  also 
observe  that  these  functions  can  be  defined  as  analytic 
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and  single  valued  In  a  neighborhood  of  real  axis  cut 
along  the  negative  real  axis  from  -a  to  0  and  along 
the  positive  axis  from  a  to  oo  .   To  start  with,  we 


•J 

suppose  U  /gh  >  1, 


I  ,  =  11m 

at 


Is^uTTe^^^ 


cosh  s  (y+h) 


•If^t   . 
e       e 


If 


^■^^-a  2  /pT  cosh  sh  /s  (a^-s^)  tanh  sh 


f , 


■  + 


The  path  of  Integration  Is  along  the  real  axis. 
The  singularities  of  the  Integrand  are  at  s  =  +  a. 

The  double  zero  at  s  =  0  of  the  denominator  of  the 

2 
Integrand  Is  neutralized  by  the  factor  s   In  the 

numerator.   We  write 

a 


I  ,=  lim 
at  , 

t->-OD 


+ 


Is^ulTe^^^  cosh  s(y+h) 


J        J  U        \  ,_  /       p   p 

-a  -r-,   r-,  J  2  /pT  cosh  sh  /s(a  -s  )  tanh  sh 


■Itf 


+ 


•Itf  \ 


ds 


11m  (l   +  I^  +  I  \ 
:-^  oo  ^  -^    ^1    a  ^ 


ds 


where  r.  Is  so  small  that  f  (s)  do  not  have  any  zero 
other  than  possibly  s  =  0  In  [-r-,,r-.].  This  Implies 
that  f  ,(s)  do  not  have  any  stationary  point  In  [-r-.,r^] 


when  U   >   gh.   Hence  I 


0(l/t)  ^  0  as  t 


oo 
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-V. 


I        =      11m 
-a 


oo  n-  -^     Isx      -Itf , 

>        b    (s+a)  e  e  ""■  ds 

t-Kc     ^         n=0     " 
-a 


"^1  1 

+     lim      /     >       c    (s+a)  e  e  ds 

t^^     '-'       n=0     " 
-a 


1-    -r. 


lim 

t-HX) 


bQ        Isx     -Itf^ 

e  e  ds 


-a 


Vs+a 


-^H 


+ 


Cq       isx     -ltf_ 
e  e  ds 


-a 


ys+a 


1-    -r. 


+     lim 
t->-oo 


r     ,oo 

/       ZZ  b    (s+a) 

-a 


n-  -^     Isx     -Itf^ 

e  e  ds 


-Pn 


o  00  ^~  ~P        ^^^        -ltf_ 

(13)  +     /  ZZ  c    (s+a)        ^   e  e  ds 

^  n=l 
-a 


-r  1 

lim      /     y       (s+a)  e  e  ds 

t->oo    ^      n=l 
-a 


The  series,   being   Taylor's    expansion   Is   uniformly 
convergent.      Hence 
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~^1  1 

lim   /  ^   (s+a)     e    e      ds 

t  ->-oo  ^   n=l 
-a 


CO    p  ■*■     n-  -^  Isx  -ltf+ 

e    e      ds 


=  11m  YH     /   Cs+a) 

t  ->-  00  n=l  ^ 

-a 


(14) 


00 


n=l 


"1        1 

lim   /   (s+a)     e    e      ds 

t  ->oo  '-' 
-a 


To  consider  these  Integrals  and,  analogously,  those  in  the 
second  parentheses  In  (15)  we  have  to  find  the 
stationary  points  of  f  (s). 


f _^  =  sU  +  ys(a2-s^)tanh  hs  ^ 


^ 


2  ^_2 


f .  =  U 


^  /T  J  (a  -3s  )  ytanh  hs  ^  h  sech  hs  ^^2   ^J> 


2  /sa^-s^ 


ytanh  hs 


for  s  7^  +  a . 

Clearly,  s  =  0  is  not  a  zero  of  f  (s).   On  the 
other  hand,  f,(s)  =  0  for  s  satisfying 


U 


/tanh  hs(sa^-s^)+  ^7"p  i  (a^-3s^)tanh  hs 


+  h  sech^  hs  (sa^-s^) >  =  0  . 
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We  write 


(15)  y  -  u/tanh  hs  (sa^-s^) 


(16)    y  =  -  1^  {(a^- 3s^)tanh  hs  +  h(sa^-  s^)sech^hs} 

and  look  for  the  intersection  of  these  two  curves. 
Equation  (I5)  has  to  be  understood  as  a  positive  branch  of 

2    2  P    "^ 

y  =  U  tanh  hs  (sa  -  s-^)  . 

This  curve  is  symmetrical  about  the  y-axis  and  s-axis. 
Also,  (0,0)  lies  on  it. 

2  ^  =  U  {h  sech^hs(sa^-s^) +  (a^-3s^)tanh  hsl 
ds  J  izz: 

ytanh  hs  (sa  -s^) 

dv 

•^  =  00   for  s  =  +  a  only 

-||  =  0   for  h  sech^hs(sa^-s^)+(a^-3s^)tanh  hs  =  0  , 


i.e.  for  s  In  [ — ,  a] .   The  curve  (I5)  therefore  is 
represented  by  the  solid  line  below. 
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We  now  consider  the  curve: 


y  =  -  "I  /^  [  (a^-  5s^)tanh  hs  +  h  sech^hs(sa^-s^) } 


y(0)  =0  and  y  <  0  in   (0,  -^) 

/5 


-  a   tanh  sh  >  0  ; 


we  know  that 

y(a)  =  0  for  (I5) • 

To  dlsglntuish  the  two  curves  we  wrlte^ 


y^  =  U  /tanh  hs  (sa^-  s^) 

yg  =  -  -|/^  {  (a^-  3s^)tanh  hs  +  h  sech^hs(sa^-s^)} ; 


then 


Thus 


ds         -y  T   y^    -^2 


-^   =  0  implies  yg  =  0 


Since  dy^/ds  =  0  for  s  <  a,  it  follows  that  the 

curve  (16)  must  enter  the  loop  of  (I5) .   Also  since  we 

know  that  yp(a)  >  0,  (I5)  and  (16)  must  intersect. 

Hence  f,(s)  must  have  at  least  one  zero.   Because  of 

f,(s)  being  analytic  in  the  open  interval  (-a,a),  there 
+ 
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can't  be  more  than  a  finite  number  of  zeros  for  f  (s). 
On  -a  _<  s  j<  0,  f_j_(s)  has  a  root  only  if  U  <  gh  and 
in  that  case  it  is  in  (-a//3,  0). 

Since  f^(s)  =  -  f_(-s)  =  0,  f_(s)  has  the  same 
number  of  stationary  points  situated  symmetrical 
to  those  of  f  (s)  on  the  negative  real  axis. 


In  order  to  consider  the  integral 


+ 


ds,  we  break  the  interval  in 


-a 


such  a  way  that  the  first  subinterval,  say,  (-a, 5)  does 
not  have  any  stationary  point  of  f  (s)  in  it,  and  each  of 
the  subsequent  ones  has  only  one  such  stationary  point 
in  its  interior.   Let  one  such  interval  be  (a,p)  such 
that  a   <  7   <  ^   and  f ,(7)  =  0. 
If  m  is  its  multiplicity. 


n-  -p   isx  -itf^ 
(s+a)     e    e      ds 


a 


(7+a) 


n — ^  17X 


2ir 


t|f^+^(7) 


m+1 


+  0  t 


/  -  TnTT  \ 


0  as   t 


00 


Similarly,  the  contribution  from  other  stationary  points 

to  the  first  term  of  the  asymptotic  expansion  of  this 

integral  vanishes.   Finally  to  consider 
5        1 


(s+a) 


n-  -^  isx   -itf-i- 


ds .   Let  s+a 


-a 


Then 
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/5+a  2 

/A  e  e  exp    l-ltiU{A+a    )+ 

0 

+  A  /( A^-a )  ( 2a-?v^ )  tanh  h ( A^-a    }    2A  dA    . 

The  major  contribution  to  this  nttegral  comes  from 
the  stationary  points  of  U(a+A^)  + 


/     2         2  2 

+  A  /{a  -a)(2a-A  )tanh  h(A  -a)    I.e.  from  the  values 

I    O 

of  A  satisfying  f  (A  -a)  =0.  As  discussed  before, 

2 
this  equation  does  not  have  any  root  because  A  <  a 

2 
and  U  >  gh.   Hence  the  first  term  In  the  esymptotlc 

expansion  of  this  integral  equals  0(l/t)  ->■  0  as  t  ->-oo. 

Similarly,  the  Integrals  in  the  second  parentheses 

in  equation  (13)  approach  0  as  t  ->  co  .   Finally,  we  have; 

b^    isx  -itf+ 

e    e      ds 


lim 

t— >oo 


>/s+a 


-a 
-r 

-L 

+   /   e    e      ds 

-a 


^        Cq    isx  -itf4. 
/s+a 


2 

If  we  carry  out  the  substitution  A  =  s+a,  the 

resulting  integrals  are  of  the  same  type  as  discussed 

earlier.   We  can,  therefore,  conclude  that  they  vanish 

for  t  ->■  CO  .   Hence  I^->-Oast->-oo.   Similarly, 

—a 

I  -*■    0  as  t  -»-  00  . 
a 
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OD 


"00 


llm 


Is  UTT  cosh  s(y+h)    e 
p    cosh  sh 


isx 


s  U^+s(-^  s^-g)tanh   sh        2J-  s(s^-a^)tanh   sh 


•tf 


(1 


+ 


.(i: 


•tf 


(1) 


f 


(1)  J 


ds    , 


where 


f i""- '    =   1  s  U  +     /s(s^  -  -g)tanh   sh 


Let 


(18) 


00 


oo 


Is^uTTcosh   s(y+h)e^^^" 


a        2  pcosh   sh  /s  — (s    -a    ) tanh   sh 


■tf^l^ 

TTTT 


ds 


For  the  convergence  of  the  Integral,  we  observe 
that  the  singularity  of  the  integrand  at  the  lower  end 
point  is  an  integrable  singularity  and  so  is  inconse- 
quential.  For  convergence  at  the  upper  end  point 
has  been  discussed  before  i.e.  it  is  made  possible 


by  the  presence  of  the  factor  e 

"^  2 
Let  s  =  a  +  t-'z  .   Then 


-s' 


m 


17(3) 
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00 
00 


.5„2,2 


.3„2 


3   2. 


lU(a+t-^z^)     TT(a+t^z    )cosh{(a+t^z    )(y+h)]e 


lx(a+t  z    ) 


b    _2p    cosh(a+t^z^)h[-(a+t^z^)t^z^(2a+t^z^)tanh  h(a+t\^) } -"-/^^ 


-t{lU(a+t^z^-  /(a+t^z^)-t\^(2a+t^z^)tanh  h(a+t^z^)}" 
o  P 


iU(a+t  ^^)-  /(a+t  ^^)-  t-^z^(2a+t^z''^)tanh   h(a+t^z^) 


2t  z  dz 


-lUta  |rP  -iUt^z^+ix{a+t^z^) 

i  e  t    /   Mz.t)  .  e  dz 

0 


where 


<l)(z,t)  = 


.5.2^2   ^,^^5  2 


U(a+t^z    )^  gCa+t-^z    )    cosh   [  (a+t^z    )  (y+h) } 


3.2, 


p    cosh(a+t\^)h{-(a+t\^)    (2a+t\^)tanh  h(a+t^z^)j 


V2 


1 


iUCa+t^z^)-  /-  t5z^(a+t^z^)(2a+t5z^)tanh  h(a+t^z^) 


J 


f-Ca+t^z^)^   t^^/^z/-(z2+  a)(^2^  ^a^^^^^  h(a+t^z'^) 
\e  -e  P  t  ^ 


^  J 


(1>1  •  <t>2  '  ^3 
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The   function   (t)(2,t)    is   bounded    for  all  positive   t   and 
is   differentlable   in   the   interval    [0,oo). 


^2 


lU(a+t^z^)-  Z-^"  t^z^(a+t-z^)(2a+t-z^)tanh  h(a+t\^) 

-iUCa+t^z^)-   /|  t5z^(a+tV^)(2a+t^z''^)tanh  h(a+t\^) 
ti . 


u2(a+t^z^)2   +^  t^z^(a+t^z^)(2a+t^z^)tanh  hfa+t^z^) 

. ,       .„  3/2   -iUta 

=      -   i   (tig  -  (^2    •    it        e 

ho  "^   P 

/i(-Ut   z   +xa+xt^z    )  ^/p     -lUta 

<t)(z,t)    e  dz    =   t^/      e 

o  i(-Ut   z^+xa+xt^z^)        3/2-lUta/T,      i(-Ut   z^+xa+xt\^) 

/   (P'    e  dz-lt      e  /i"e  dz 


where 

We  have,  therefore,  to  consider  the  two  integrals 

-lUta  ^^2^/<t>'|  K-Ut^z^+xa+xt^z^) 
t   ^  <    ,e 


e    L '  ■   ;  ^    ;  e  d  z 


0 
iUta  3/2^7^/(1)'")  l(-Ut^z^+xa+xt-z^) 
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/t7p 


where 

- lUt a J57p  (■(})•  ]  l(-Ut^z^+ax+xt^z^) 
J,  =   e    t^^/   ^    >  e  dz  . 

The  sign  of  (|) ' ,  cj)"  depends  upon  the  sign  of  the  function 
we  write  J^^e-^t-t'/^y ''^  *■  ^K-Ut^z^+ax+xt 'z^)^^  _ 

o 
The  function  cj) '  being  analytic  will  have  isolated 


and,  therefore  finite  In  number,  zeros  In  the  Interval 
[0,  /T/p] .   Also  the  number  of  points  In  [0,  yT/p  ]  at 
which  <j) '  attains  maxima  or  minima  will  be  finite.   These 
two  observations  together  imply  that  it  is  possible  to 
divide  [0,  )/T/p  ]  into  a  finite  number  of  such  Intervals 
such  that  (j) '  is  monotonlc  in  each  of  them  and  maintains 
the  sign.   Let  one  such  Interval  be  (a,p)-   The  contri- 
butlon   to  J,  from  this  interval  is: 

-iUta      lax    ^/o /?  l(-Ut^  +  xt^)z^ 

I  =   e  etM'e  dz. 

a 
There  are  four  cases: 

(1)  <t) '  is  positive  and  monotonlcally  decreasing  in  (a,^) 

(2)  ({) '  is  negative  and  monotonlcally  increasing  in  (a,p) 
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(5)   "f) '  IS  positive  and  monotonlcally  Increasing 
(4)   (f) '  Is  negative  and  monotonlcally  decreasing, 


Case  1 


-lUta  lax   3/2, 


I=e  e        t/(j)'e 

a 


l(-UtVxt^)z^ 


dz 


lax-lUta  -^/rp 
+  e  t^^    (t)'{a)/    e  dz    , 


^    l(-UtVxt^)z^ 


a 


for 


a 


and 


l(-UtVxt-^)^^ 


/?    l(-UtVxt^)z^ 


dz  = 


^    <  P 


1(-Ut>xt")a 


2U(-UtVxt^)  2ia(-UtVxt^) 


+  0( 


(-Ut    4-xt  0 


)  for  a   ?^  0 


a 


^(Ut'^-xt^)  I 


^1/2        1       . 


+  0(  r^ — — 

(Ut   -xt") 


for  a   =  0, 


Also   llm  <}) '  (a  ,  t)  =  finite.   Hence  I->-0  as  t->-ao. 
t^-00 

In  case  (2),  we  consider  -  I  =  gi(ax-Uta ) ^3/2 


(-<i)')e" 


l(-UtVxt")z^ 


dz  .   Since  -([) '  Is  positive 


-  a 


and  monotonlcally  decreasing.  It  follows  from  case  (1), 
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that  -I  ->-  0  as  t  ->■  oo  . 
In  case  {^) ,    we  have. 


lax-lUta/?,  ^   i(-UtVxt^)z^ 


-I  =  t^/^e        /  f(l)'(3)-(t''(z)]  e 


dz 
a 

y2-l(ax-Uta)  P  i(-Ut^+xt^)z^ 

■t   e      ^   /  <l>'(p)  e  dz  . 

a 

In  the  limit  t  ->-  cx) ,  the  first  integral  on  the 
right  vanishes  as  in  case  (1) . 

From  the  asymptotic  approximation  of  the  second 
integral,  it  is  clear  that  the  second  Integral  also 
tends  to  zero  as  t  ->-  oo  .   Case  (4)  reduces  to  case  (3) 
if  we  consider  it  as 

iax-iUta  ^/P /?     l(-Ut^+xt5)z^ 
-I  =  e        t'^/'y(-(t,')e  dz 

a 

and  therefore  can  be  treated  likewise. 

As  (a,p)  was  any  arbitrary  subinterval,  the  same 
can  be  concluded  for  all  such  Intervals.   Hence 
J^  ->  0  as  t  ->■  00  .   Similarly  J^  ->■  0  as  t  ->■  oo  . 

We  now  consider 

-lUta  /  p°  f  (k  t  7  i(-Ut^z^+ax+xt^z^) 

VWp 

/  iax-iUta /7         i(-Ut+xt^)z^ 
=  It-'/'^e        /  (t)(z,t)  e  dz  . 

/T/p 
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00 


Jgl  ^  t^^^  f  \H^,t)\  dz 


/tTp 


(1)1   <  1^1 


.3.2, 


^     =  Ut^(   a   ^  ^2)2        g   cosh[(a+t^z    Hy+h)} 
^  t^  p    cosh    (a+t^z    )h 


5   2,, 1/2 


t^[-(^  +z^)(^  +z2)tanh  hCa+t-'z^)] 


For   large   t. 


p  (a+t-'z    )y  T     /-  ^ 


p  '   T  —  p  y 


T 


where    |g|   _^  M/U   for  all   z   In   [0,oo). 
-t^z^      VW-9    z^t^l/^ 


(t)^|    %  e 


oo 


Jo        < 


^^^'        /'t^z^MJ^e 


,-,      -t^zV/T7F   z^tll/2 


■2'    -       p  J     "    "    "^T 

yr/p 


dz 


.9/2 


oo 


-  ^  M  ^V^  ^ 


-  A^^^-  /t7^  A^ 


dz    . 


/tTp 

Let  z  =  /tTp  -2  ,   dz  =  /tTp  -^ 
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^'^'^    (^)    llm      r     i-    e 


2      t^(T/p)2(i^-i^) 


P  P 


I>>oo 


M 


T   t 


y2 


6,„  ,  .2,«3  A 


llm      /  {{t^{l)^{^Jl^-k£^)    e 


t:^ 


6£ 


Tt 


M 
57^ 


llm 
L->-oo 


-   k    f    -^ 5 6£ 


M 


llm  U-  u(   ^^  +  1)>  , 


-T2      -J-iin  AX-  m    -^_j.'-j-    +  1 ;  |>  ,      where  0   <    n  <   1 

1  U  J_r"^KX)  I  j 


1-a 


0        as   t  ->►  00 


T   t 


I.e.    I„    -»-  0  as    t  -*■  00  .      We   then   consider 


"00 


00 


00 


3%"TTcosh  s  (y+h) 


-tf 


(1) 


+ 


a      2p    cosh   sh  /s  — (s    -a    )tanh  sh 


TTT 


+ 


^^ 


00 


I   < 

OD  '  — 


s^U 


cosh  s(.v+h) 


a   2  /pT  cosh  sh  ys(s  -a  )  tanh  sh 


/   2  T 
-t  /s(s  —  -g)tanh  sh 


f 


(1) 


ds 


This  integral  Is  unlformally  convergent  for  all  t  >  0. 
Hence 


00 


11m 


sU ITT  I cosh  s(y+h) 


^"*'®  a   2  y^T  cosh  sh  /s(s^-a^)  tanh  sh 


00 


a 


lim 
t^>oo 


2  T 
-t  /s(s  —  -g)tanh  sh 


f 


(1) 


ds 


-J     ■ — - 

s  u|TTIcosh  s(y+h)  e 


-t  /s(s  —  -g)tanh  sh 


2  /pT  cosh  sh  ys(s^-a^)tanh  sh  [f^"'"^  | 


ds 


0 


Hence  I   ^  0  as  t  ->■  oo  .   Similarly  we  can  prove 


that  the  time-dependent  part  In  I 


-GO 


0  as  t  ^^  oo  . 


It  follows  from  (11),  therefore,  that  time-dependent 

terms  approach  zero  In  the  limit  as  t  -*- oo  ,  I.e.  In 

2 
case  U  >  gh,  the  steady  flow  exists.   In  order  to  find 

the  relevant  radiation  conditions,  we  collect  the 

remaining  terms.   We  have  then. 
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00 


Is  uTTe    cosh  s(y+h)  ds 

^'^     -00   p  cosh  sh(s^U^-/-  s(a^-s^)tanh  sh) 

r 


We  have  to  find   limit  ^    .      Consider 

x->  +00   ^ 
-00 


00 


(19) 


lU  cosh  s(y+h) 
cosh  sh 


CO 


-00 


1  ^^r    iU  cosh  s(y+h) 
J   J  cosh  sh 

0 


cosh  s(y+h) 
cosh  sh 


eS(y4-h)^  3-s(y+h)   ^  ^3y  ^^^-2s{y+h) 


^sh  ,  ^-sh 
e   +  e 


1  +  e 


-sh 


We  observe  that  s  >  0  and  y+h  >  0 


sy  i+e-2s(y+h)     3y 

<  2e 


1  +  e 


■2sh 


since  y  Is  negative,  (19)  converges  absolutely. 
The  remaining  function  is 


(20) 


Tr(s)s^ 


2  2         2  T 
s  U   -  s(g-s  — )  tanh  hs 
F 


in   (-00  ,00  ) 


The  singularity  of  the  function  is  a  double  singularity 

2 
at  s  =  0  -which  is  neutralized  by  the  factor  s   in  the 


5^ 


numerator.   This  Implies  that  the  function  (I9)  Is  bounded 
Hence 

r   Is  U  n  cosh  s (y+h) 

-00    p  cosh  sh(s  U  -  s(g-s  -)  tanh  sh) 

is  absolutely  convergent.   By  the  Rlemann-Lebesgue  lemma, 
therefore. 


00 


Isx 


1^        [      i3  uTTcosh  s(y^h)  e- ^^  ^  ^ 

^l"*^_cn   P  cosh  sh(s  U  -  s(g-s  -)  tanh  sh) 


Also  If  we  had  worked  with  i   instead  of  working  with  i  , 

y  "=      ^x' 

the  only  difference  is  that  we  would  have  sinh  s(y+h) 

instead  of  cosh  s(y+h).   Exactly  the  same  analysis  holds 

in  this  case,  which  implies  that   lim  \>      =  0,  i.e.  the 

|xI-»-oo  ^ 
disturbance  dies  out  at  +  00 .   Hence  the  conclusion  is 

2 
that  in  case  U  >  gh,  the  steady  solution  does  exist  and 

the  corresponding  radiation  conditions  are  that  the 

disturbance  dies  out  both  upstream  and  downstream. 
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Case  II    U  =  gh 
Again 


r  a 


<l>(x,y;t) 


/2? 


A{s,t)  cosh  s(y+h)  e^^^   d 


-a  OD 

+  I     A(s,t)  cosh  s(y+h)  e^^^  ds  +  rA(s,t)cosh  s(y+h)s^^^ds 
-co  a 


/27r 


a 


r   IsUTTcosh  s(y+h)e 
p  cosh  sh 


isx 


-a 


_s  U  -s(g-s  — )tanh  sh 


+ 


-Itf^     -ltf_ 
e         e 


/  2  T 

2  y(g-s  — )s  tanh  sh 


f 


+ 


f 


ds 


+ 


-a    cx)' 


+ 


CO   a 


■sUTTcosh  s(y+h)e 


ISX 


p  cosh  sh 


.s  U  +s(-g+s  — )tanh  sh 

r 


e  e 


2  /s(s^  -  -g)tanh  sh 


f 


w         ;rTT 


f 


+ 


'J  J 


ds 


f ,  =  sU  +  ys(g-s^  -)tanh  sh 


(1)  / p-pp 

f^   =  IsU  +  /s(-g+s  -)tanh  sh 

r 


+ 


^^a  ^at    ^^00  ^OD  t   ^^-00  ^-cDt 
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We  write 


a 


,t  ^      U  r      TTscosh   s(y+h)    e^"^        ,    ""f^     -Itf 

^^  /2Tr     _-'        /pT     /s(a   -s    )tanh   sh 


-a 

The  main  contribution  to  the  asymptotic  expansion 
of  the  integral  cj)  .    comes  from  the  stationary  points  of 
f  (s)  and  from  the  endpoints.   We  have  already  proved 
in  a  similar  integral  that  the  contribution  from  the 
endpoints  tends  to  zero  as  t  -*-  00  .   For  the  contribution 
from  the  stationary  points,  we  study 

/  =  -^     f      Tts  cosh  .  (y-H.)  e"'"   /  ""^  -i^^'^as 
1^"^^     -a   /pT  ys(a  -s  )  tanh  hs  ^ 


f_(s)  has  a  double  zero  at  the  origin  and  f^(s) 

does  not  have  any  zero  at  the  origin.   Let  us  suppose 

that  all  other  zeros  of  f+(s)  are  simple.   Then 


t 


+   2U 


at 


/27rpT    La  -/h 


— 2Z )l/5^  0(1) 

t|f'"(0)l         t^/^ 


i.e . 


4    «"''"> 

(})      0(t   ■^)    which  ->  00    as  t  -*  00  . 
a  T3 


I 


If  f  (s)  has  a  zero  of  order  higher  than  2.  the  rate  at 

which  (|)  ,  -*- CO  will  be  accordingly  higher, 
at; 


39 


On  the  other  hand,  the  contribution  from  I    ,,  I^  ^. 
has  already  been  seen  to  tend  to  zero  for  t  ->  oo  . 
Hence  in  the  critical  case,  there  does  not  exist  a 
steady  motion.   The  velocity  potential  ->■  oo  at  least 
as  t  ^  ^.  For  any  discussion  of  this  case,  we 
have  to  revert  to  linear  theory. 


ho 


Case  TIT:   U^  <  gh 


Isx 


^TT   =  —-Z     /  A(s,t)  s  cosh  s(y+h)e    ds 


*    -oo 

-a   a   00 
ff  /  +  /  +  /  ^  A(s,t)s  cosh  s(y+h)e    ds.] 


2Tr 

-oo   -a   a 


~   'I-c  -  ^a  -  Ice  > 


p  Isx 

Ig.  "   /  A(s,t)s  cosh  s(y+h)e    ds 


-a 

a 


i  s^  U" 


[ 


Ji    p  cosh  s  h   s^ir-s(g-s^  -)tanh  sh 

a.  p 


+ 


/  -if ,t       -If  t 
^  fe   +    _   e   - 

2js(g-s2(T/p))tanh  sh    V    ^+  f_ 


isx 
'Cosh  s(y+h)e    ds 


f ,  =  s  U  +  /s(a^-s^)tanh  sh  /^ 
±  -  ^  p 


Tn  this  case,  we  know  that  both  f_^  has  a  simple  zero  each 

at  s  -  0.   Besides  f_(p)  =  0  =  f  (-p). 

We  write  I^  =  i   +  I  ^  , 
a    as    at  ' 


where 


a 


2,,  ^    „lsx   ,   ,     .      -if+t    -If  t 


I   =   r   i  s ^  TT  e   cosh  s(y+h)     /O   "^    e     ,^ 
at    J      2p  cosh  sh   ^—  ^ f f 'ds 

-a  /s(g-s  -)tanh  sh       ^       " 


hi 


2 
Owing  to  the  presence  of  the  factor  s   with  numerator, 

s  =  0  is  not  a  singularity  of  the  Integrand. 

In  order  to  consider  the  singularities  at  s  =  ±  P,  we 

deform  the  path  of  Integration  at  s  =  p^-p  with  two 

semicircles  in  the  lower  half  plane  as  shown  in  the 

figure . 


-P      ,       P 


r 


>-  I  V  [      , V 

0 


f  _  =  sU  -  /(gs-s^  -)  tanh  hs 


2  T, 


r 

-1    (g-3s  — )tanh  sh+h  sech  hsCgs-s-^— ) 
=  u  -  2  [  ^ 7z:=^ ^  1  • 

/(gs-s^  — )    tanh   sh 


Also 


PU^  =    (g  -  P^  ^)    tanh  ph 
,  T      pU^-2p^  -   tanh  ph+h  sech^ph  p^U^coth  ph 

f  _(P)  =  u  -  I  ( e ^ } 

-  ^  {U2   +  2p  ^   tanh  ph  -  -i^^} 


W  fU^f^   -   slnhlhp^    -^  2p  I  tanh  hp}      >   0 
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Similarly 


f   (s)    =   sU  +    /(gs-s^  -)    tanh  hs 
,  n       {pu^-2p^  ^   tanh  hp    +  h   sech^hp   p^U^coth  ph} 

^f-P)=  u  -  A  e p^j 

which,    before,  is  >  0 . 

Over  the  semicircles  s-p  =r-,e   ,  s  +  p  -   re  ^, 
r,r^  small  and  tt  <  0,   cf)  <  27r,   f_(s)  and  f  (s)  behave 

like  f_|_(p)(s-p)  and  f|(-p)(s4^)  i.e.  like  fjp)re^®, 

'      id) 
f  (-p)r,e  ^  respectively. 

That  is,  the  real  parts  of[-if_(s)},  {-if_|_(s)1  are 

negative . 

We  now  have: 


I      r  1  s  UTT   e^^^cosh  s(y+h) ^e ^-—-"ids 

at    J   p  cosh  sh   ^— p-pp f ,       f_ 

r  2/s(g-s  -)  tanh  hs 


■itf^ 


r  i  s^  uTT   e^^^cosh  s(y+h) 

J      p  cosh  sh    / p-fp f 

r  2/s(g-s  -)tanh  hs 


ds 


2  „=F\ 


i  s  U  I    e   cosh  sfy+l      e_ 


-itf 


ds 


J      p   cosh  sh   f- p— Tp f_ 

r  2js{g-s     -)tanh  sh 

We  consider  the  first  integral  over  the  semicircle 

s  +  P  =  r-.  e  ^  and  the  second  integral  over  the 

iO 
semicircle  s-p  =  re   .   The  two  integrands  do  not  have 

any  singularity  over  the  semicircles.   They  are,  in  fact, 


^5 


analytical  function  there  at  and  so  bounded.   Hence 

P    -     .  -Itf 

"^  1  s^  U  TT   e^^^'cosh  sfv+h)     - 


e 


J       id)  P  "^"^^  ^^      r 2~T 

!+p=r,e  ^  2/s(g-s  — )tanh  s 


h 


^ 


ds 


-Re(tf ,) 
M  C      ^   Ids 


s+p=r^e 


0  as  t  ->  00 


where  M  is  a  bound  of  the  coefficient  of  e~"^^-^+  in  the 
Integrand  on  s+p  =  r,e  ^ ,    ir  <   ^   <   2tt . 

Similarly  the  second  integral  also  vanishes . 
On  the  straight  parts  of  r  and  also  for  I_|_^  ,  we  can 
proceed  as  before  in  the  first  two  cases,  to  prove  that 
as  t  ->  oo  these  integrals  vanish. 

Hence  the  steady  motion  in  this  case  does  exist. 

To  find  the  radiation  conditions  at  +oo,  we  consider 
the  remaining  integral, 

^  _   1   r   i  s^  uTT    cosh  s(y+h)e^^^     , 

^  ^  ,/2?  J   P  ^°^^  ^h   s2u^-s(g-s^  ■^)tanh  sh    * 
'   -oo  p 

Let  us  suppose  that  x  <  0.   The  zeros  of  the  denominator 
are  at  s  =  0  (double  zero),  and  s  =  P , -p •   Owing  to  the 
presence  of  s   in  the  numerator,  the  only  singularities 


rom 


are  at  s  =  +  p .   We  again  consider  the  contour  p  f 

-00  to  00 .   On  this  contour  the  coefficient  of  e    is  again 


41^ 


analytic  function  and  so  bounded.   On  the  two  semicircles, 

jLg^    ixr^(cos  (|)+1  sin  ^) 
Q  =   e  7r<(j)<27r 

and 

Isx    ^-^^    ^'^°^  ®"^^  sin  0) 
®"^  7r<e<27r 

This  =>  that  for  x<  0,  the  real  part  of  the  exponent 

on  these  circles  Is  negative.   We  can  then  argue  as 

before  to  conclude  that  for  x  ->oo,   I  vanishes 

s 

On  the  straight  parts,  the  Integrand  has  no  singularities 
and  the  Integral 


o     - — 

(ni  )  r  i  s  U  TT  cosh  s(y+h)   r  1  ,  _^ 

^^^W  p  cosh  sh  L  -2-2   ;    ^^ J  ds 

s  U  -s(g-s^  -^)tanh  sh 

r 

Is  absolutely  convergent. 

Hence  an  application  of  Rlemann  Lebesgue  lemma 
shows  that  the  above  Integral  I  vanishes. 

It  follows,  therefore,  that  the  disturbance  vanishes 
upstream. 

Finally,  we  consider  x  >  0 
Ks)  = 
(22)    ^   r?         Is^  U  TTcosh  s(y+h)e^^^ 


2tt2   /    2  T, 


V2Tr_^         p  cosh  sh[s  U  -s(g-s  -)tanh  sh] 

P 

The  Integrand  now  has  poles  at  s  =  p,  -p. 


ds  . 
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We  deform  the  path  of  integration  at  s  =  p,-3 

by  drawing  semicircles  In  the  lower  half  plane   On 

the  straight  parts  of  the  contour,  argument  similar 

to  that  given  above  shows  that  l(s)  =  0  for  large  x. 

the  only  contribution  that  we  get  is  from  the 

semicircular  parts,  say,  B  ,  B_ .   To  get  that  we  first 

complete  the  two  circles  so  that 

.^^    isr  (cos  <|)+i  sin  i?) 
e    =  e  on  B_ 

isr  (cos  0+  i  sin  0) 
=  e  on  B   . 

On  0  <  0,  (j)  <  TT,  the  real  part  of  the  indices  of  the 
exponential  is  negative.   Hence  when  x  -»■  oo  ,  I(s)  ->- 0 . 
This  shows  that  the  integral  over  the  semicircle  is 
the  same  thing  as  that  over  the  whole  circles. 

On  the  other  hand,  since  the  integrand  is  analytic 
function  except  for  two  poles  at  the  center  we  can  use 
Cauchy's  residue  theorem  whence  It  follows  that: 

ls^U(lT(  s  )  e^^P+tt(  -s  )  e"^P  )  (  s-p  )  cosh  s  fer+h) 

I  *=   -i.V^  lini   { 2~2 2~T ^ 

s->p     cosh  sh(s  u  -s(g-s  — )tanh  sh 

where 


H  =  11m 


Jj:g_ 


2.  2       2  T 
s->p   snj  -s(s-s  — )tanh  sh 
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p  Is  given  by  the  equation 

Pl^^  =  (g-p^  ^)    tanh  ph  . 
For  large  x,  therefore, 

which  obviously  does  not  tend  to  zero  for  x  ^-  oo  . 
The  disturbance  therefore  does  not  die  downstream. 

We  can  evaluate  H  explicitly. 


H  =  11m 


_§j:e_ 


?  2       2  T 
s->p  s  U  -s(g-s  — )tanh  sh 


2pU^-(g-3P^  -)tanh  ph-Cgp-p^  ^)h   sech^  hp 

^  1 

2(g-p^  -)tanh  ph-(g-3p^  -)tanh  ph-(gp-p^  -^)h  sech^hp 
P  P  P 

^ 1 

(g+p^  -)tanh  hp)-(gp-p^  -)h  sech^l^ 

P  r 

Hence,  when  T  -»►  oo  ,  p  ^  0;  H  varies  as  1/5  and  (|)^  ^  0 . 
The  disturbance,  therefore,  though  not  dying  downstream 
In  general  for  any  finite  T,  -»-0  when  T  ^  cd  . 
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We  can  now  summarize  the  results.   First, 

steady  state  exists  both  for  U  /gh  <  1.   For  the 

2 
critical  case  i.e.  U  =  gh,  there  is  no  such  solution. 

Indeed,  the  velocity  potential  grows  with  t  at  least 
as  t    and  becomes  infinite  as  t  ->■  oo  .   In  case  of 
supercritical  case,  the  radiation  conditions  are  that 
the  disturbance  should  die  both  upstream  and  downstream. 
In  the  subcrltical  case,  the  disturbance  tends  to  zero 
upstream,  but  does  not,  in  general,  do  so  downstream. 
When  it  does  not  tend  to  zero,  the  velocity  at  infinity 
is  given  by  (23) .   If  we  take  T  =  0,  we  get  the  same 
result  as  obtained  by  Stoker  [9]-   On  the  other  hand, 
the  dependence  of  (j)   on  surface  tension  T  is  inversely 
proportional.   When  T  -*- 00  ,  the  disturbance  dies  down- 
stream as  well . 

Since  gravity  being  a  finite  force  is  negligible 
in  comparison  with  arbitrarily  large  T,  the  above 
result  holds  even  for  non-gravitating  fluids. 


48 


Chapter  II 
The  Solitary  and  Cnoldal  Waves 

Formulation  of  the  Problem 

We  consider  two  dimensional  flow  of  liquid  in  an 
Infinite  channel  of  finite  depth  d.   The  assumptions  on  the 
liquid  and  the  flow  are  the  same  as  those  in  Chapter  I. 
The  pressure  on  the  free  surface,  however,  is  supposed  to 
be  Pq(x) . 

Let  u,  V  be  the  horizontal  and  vertical  components  of 
velocity  at  a  point  (x,y)  at  a  time  t  after  the  disturbance 
has  been  generated.   The  equations  governing  the  flow  are: 


U,   +UU   +VU   =-—  p 

t      X      y     p  ^ 


V.  +UV   +VV   =-—  p   -g 
t      X      y     p  ^y   ^ 


V   =  u 
X    y 


The  equation  of  continuity  gives 


u  +  V   =0 
X    y 


The  kinematlcal  conditions  are 


T]^  +  u  T]^  =  V     at  y  =  Ti(x,t) 


V   =  0     at  y  =  d 


^9 


The  discontinuity  at  the  free  surface  is  given  by 


r,"^ 


where  p  (x)  is  the  pressure  of  the  atmosphere  at  a  point 
(x,T^)  just  outside  the  liquid  and  p(x,rj')  at  the  corresponding 
point  inside  the  liquid.  T\{y,,t)    is  free  surface  elevation. 
We  suppose  that  P„(^)  — ^  0  as  x  — >  -  oo .   The  problem  then 
is  to  obtain  the  equation  governing  this  unsteady  wave 
motion  by  using  the  stretching  and  expansion  procedure  of 
Keller  [ 7  ]  and  Frledrichs  [  2  ] •   This  equation  is  then  to  be 
solved  for  particular  pressure  distributions  on  the  free 
surface. 

We  use  the  following  transformations  to  change  the 
independent  and  the  dependent  variables  in  the  above 
equations . 

^   =  '^  y  =  Ph     u=  \/gh  u     V  =  \/gh  1^ 

T  =  v/gh  00^  h^  t     Ti  =  hY     d  =  hH     Cf  =  co^h^ 

P  =  gph  TF       Pq  =  gph  0     TTg 


where  h  has  the  dimension  of  length  and  cd  has  the  dimension 
of  inverse  length.   The  new  quantities  a,  p,  u,  v,  tt,  tTo, 
y,  T  are  dlmenslonless .   We  have  then: 
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u^    ^  v/gh  a.   u^  \  =/eh  u^    ■    ^ 


u ,    =  /gh    •  v/gh  o)-^   h      u     =  gcD-^   h-^   u  v      =  v     i/  ^ 


oj 


y/gh  ~  2    ,  2 


Py    =    pgTTp  Tl^    =   Y^  CD^    h^  v^  ^x    "    ^   ^a 


^    ^^ -  -    ■-      O) 


T   =  gph^   Tq 


substituting  these   In  the   above   equations  we   get, 


(T  u^  +  vp   =  0 


2 

<^     u^   +    cTu  u^   +  V   Ug  +    cTtt^^  =  0 


^    v^+cruv^+vvQ  +  cr7rQ+(r=o 

T        a      p      p 


^a  -   ^p 


TT  -  d^  TTc;  =  o-^To    at  p  =  Y 

^    d+O-Y"^)^/^ 
^     a 


V  =  0  at  P  =  H 

^^  Y^  +  cTu  Y^  =  V        at  p  =  Y 
We  assume  the  following  expansions: 
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J2.  3 

u  =  u     +  iTUt    +  'J    u„   +  <y  n-,   +   . . 
o  L  d  ^ 


2  3 

o  1  2  3 

2  3 

TT    =    TT       +  CtT-,      +    <J    Tr_  +    •CT    TT^     +     .  .  . 

o  1  2  3 

Y   -   Y      +  ■-TY-,     +  CT^Y^  +  cr^Y^    +    .  .  . 
o             1                2  3 


Substituting  these  for  u,  v,  tt,  Y  In  the  above  equations 
we  get. 


O'fu   +  c'Ju,   +  f5"  u„   +  (J-^u^   +  ...  1  +  V  Q  +  crv-io  +  <y^Vr.a, 
"■  oa     la      2a      3a      -'     oP     ip      2P 

+  '^■^v^o  +  .  . .  =0 

2  2        3 

J  [u   +  aun   +  J  u^  +  cT^u^  1 
■■  OT     It      2t      3T-' 

+  ^yfu      +  cTUt     +  (T    U„    +    .  .  .  1  [u         +  'J-Ut       +  cT  u„      +  cr  u^      +.  .  .1 
'■    o  1  2  -'  '-    oa  la  2a  3a  ■' 

+   [v^+cTv-j^+o  v^+cTv-^H-.  .  .  ]  [u^g+^Tu^g+cru^o+j^u^g   +   ...] 

2 
+  (^[Tr      +0'7r-,    +^  7r„      +...1=0 
'-    oa        la  2a  -' 

+^[%^V^^S+^S   +    •••^f%a+^^la+^^^2a+''^Sa   +    •••] 
+    [v^+o"v-j_+»^V2+^  v^   +    .  .  .  ]  [v^g+Ov-j^g+fi^Vgo  +  'f-^v^g   +    ...] 

+  ''^['^^p+<J-^-^n  +  0-W       +6^TY  +     ...]     +    .J- =    0 
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"^^^^2aa  "  2  ^oa  ^laa  "^  "5~  ^oa  ^oaa  "  ^^oa  ^oaa  ^la^ 

-h^-^fy     -  —  Y    Y     +  is.  Y    V  +  is.  V     V^   V 

^  ^-"^aa   2   oa  ''2aa   ^  ""oa  naa  ^  2  ^oaa  ^oa  ^laa 

-  2^   Y     Y    -  2.  Y^   Y  -  ^Y    Y    Y    ^  +     1 

lb   oaa   oa    2   la   oaa  -^  oa   2a   oaa^  -r  .  .  .  j 


at  p  =  Y^ 


Vq  +  crv-j^  +  cr  V2  +  cr^v^  +  ...  -  0       for  p  =  H 

^(YoT-^YlT+^Y2T  +  •••)  ^'^(Yoa+'^la+^Y^.+A^^  +  ...)  ^ 
(V^^l+^lV^  ^^{u^+Y^u^P+Y^u^P+Iy^  u^pp)  + 
+  cr^(u3+Y-LU2p+Y2U^p  +  l  y2  u.pp+Y^u^p+Y^Y^u^pp 


Y? 


=  v^+,5tVi+Y^v^p)4-d2(v2+Y^Y^p+Y2V^p  +  |  Y^  v^p^) 

+^^(v3+Y^V2p+Y2V,p+  I  y2  v,pp+Y3V^p+Y^Y2V^p+  ^  v^ppp)+. 

at  p  =  Y^ 
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Equating  zero  order  terms  we  get, 

(1)  ^op  =  0 

(2)  v^  u^p  =  0 

(5)  ^o  %p  =  ° 

(^)  ^oa  =  V 

(5)  ^o  =  °  a*  P  =  Yq 

(6)  ^o  =  °  at  p  =  H 

(7)  ^o  =  °  at  p  =  Y^ 
implies  VQ(a,p,T)  =0   u^  =  u^(a,T) 

I  order  terms 

(8)  %a  +  ^ip  =  ° 

(9)  %  %a  +  ^o  ^ip  +  ^oa  =  0 

(10)  %  %a  +  ^o  ^ip  +  ^op  +  ^  =  ° 

(11)  ^oa  =  %P 
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(12)  TT,    +  Y,    7r^„    4-TY  =0  at6=Y 

1  1      op  o      oaa  au   p        Xq 

(13)  v-^   -  0  at  p  =  H 
The   first   equation  yields: 


^1   ^   -   P   %a   +  ^(«'^) 


0   =   -   H   u^^   +  f(a,T) 


(15)  Vi   =   (H   -  P)    u^^ 


The  third  equation  gives: 

^o  "  -  P  +  fi(a.T) 
0  =  -  Y^  +  fl(a,T) 


(16)  TT^  .  Y^  -  P 


Substituting   for  tt     in    (9) 


u     u        =  0 
o      oa 


%    =   ^2^^) 


We  assume  that  the  disturbance  at  a  =  -  oo  is  independent 
of  time  T. 
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Then 


u   =  const 
o 


Let 


u(-oo,t)  =  c  ,  a  constant 


c  =  u(-oo,t)  =  u   +  cTu-.    + 


{-oo,t) 


(17) 


=  c  +  (fc^    +  O      c„  + 
o     1       2 


Thus 


u  =  c 

o    o 


The  last  equation  now  yields 


(18) 


Y  =   constant 
o 


Hence  zero  order  terms  yield  uniform  flow, 


II  order  terms 


"la  +  ^2P  =  ° 


"ox  +  "o  "la  +  "l  "oa  +  ^o  "2P  +  ^1  "ip  +  ^2  "op  +  ^la  =  ^ 


^OT  +  "o  ^la  +  "l  ^oa  +  ^o  ^2P  +  ^1  ^ip  +  ^2  ^oP  +  ^ip  =  ° 


^2a  =  "2P 


We  enjoin  to  these  the  following  equations  from  the  first 
order  terms: 
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^1  +  ^1  ^op  +  '^o   ^oaa  =  0     3t  P  =  Y^ 


V2  =  0  at  p  =  H 

+  I  ^1  ^opp   at  p  =  Y^ 


Using  the  previous  results,  the  above  equations  are  reduced 
to 


(19)  u^^  +  ^2P  =  0 


(20)  u   u,    +  TT,    =  0 

o   la    la 


(21)  .^p  .  0 


(22) 


^2a  -  ^2P 


(23)  TT,  =  Y,      at  p  =  Y 

11  o 


(2^)  Vg  =  0       at  p  =  H 


(25)  u^  Y^^  ^v^  at  p  =  Y^ 


From  (21)  and  (23)  we  have  v^    =  Y^   everywhere.   Integrating 
the  equation  (20)  we  get. 
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%  ^1  +  ^1  =  ^3^"^'^^ 


(26)     i.e.    u^  u^  +  Y^  =  f  (t,P) 


In  view  of  the  assumption  that  the  disturbance  at  a  =  -oo 
Is  Independent  of  time,  f^(T,p)  In  (26)  must  be  a  constant. 
In  particular  u-.  Is  Independent  of  p. 


^2P  =   -   ^la 


^2  "  "  ^^la  "^  ^}^ict,-v) 
0  =  -  Hu^^  +  fi^(a,T) 
(27)  Vg  =  (H  -  p)  u-L^ 

From  (22)  we  then  have: 


u^„  =  (H  -  p)  u^ 
2P  laa 


Ug   =  -  I  (H  -  P)2  u^^^  +  F(a,T) 


From  (25)  and  (27)  we  get. 


u  Y.   =  (H  -  Y  )  u, 
o   la    ^      o    la 


Substituting  for  v^   in  (20) 


u   U-,   +  Y-,   =0 
o   la    la 
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The  two  together  Imply  that  either  u-,,  Y,  are  constants  or 
that 

(28)  u^  =  Y   -  H 


o    o 


and  then 


^1  =  -  %  ^1  -^  S^^^ 


Again  evaluating  this  equation  at  a  =  -  oo ,  we  find  fc(T) 
a  constant,  say  6. 

(29)  Y^  =  -  Uq  u^  +  5 

Also  since  u  =  c  ,  v-j^  =  0 

U2  =  -  I  (H  -  P)2  u^^^  +F(a,T) 

III  order  terms: 


^2a  +  ^3P  ^  ° 


^IT  +  %  ^2a  +  ^1  ^la  +  ^2  ^oa  +  ^o  ^3P  ""  ''l  ^2P  +  ^2  "ip 

+  ^3  ^op  +  ^2a  ^  ° 


V-,    +  U   V„   +  U^  V,    +  U„  V  ^  +  V^  V-,o  +  V-,  V„o  +  ^o  ^TR 

It     o   2a    1   la    2   oa     o   3p     1   2p    2   ip 


+  ^3  ''op  +  ^2P  =  ° 


3a    3P 


v^  =  0         at  p  =  H 
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^IT  +  ^oa(  )  +  ^la^^l  +  ^1  "op)  +  %  Y2a  ^  ^3  +  ^2P  ^1 


at  p  =  Y^ 

With  these  we  take  the  pressure  equation  from  the  second 
order  terms : 


^2  -  ^S  -^  ^1  ^IP  ^  ^2  ^op  +  I  Y?  ^opp  =  -  To(Yic.a  -  I  ^oa  ^oaa^ 


at  p  =  Y 
^    o 


Making  use  of  the  previous  results,  the  above  equations  can 
be  reduced  to 


(30)  U2„  +  ^3P  ^  ° 


(31)  Ut   +  u  u„   +  U-,  Ut   +  Tr„   =  0 

^■^    '  It    o  2a    1   la    2a 


(32)  u^  V2„  +  ^2P  =  ° 


(33)  v^^  =  u^p     at  p  =3  H 


(3^)  v^  =  0 


(35)  Y.   +  Y,   u,  +  u  Y_   =  v^  -  Y-,  u^    at  p  =  Y 

^•^^  It    la   1     o   2a    3    1   la     ^    o 


(36)  ^2  -  ^3  _  Y^  +  T^  Y^^^  =0     at  P  =  Y^ 
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^2P  +  "o  ^2a  =  ° 


TTgp    +    u^(H    -    P)     U^^^    =    0 


TT^    -  ^   O    -   H)2   u^^^   =   G(a,T) 


At    P   =  Y^  Tr_    =  7r„    +  Y^    -   T      Y^ 

o  2  S  2  o      laa 


^   (H   -   Y    )2  +  G(a,T) 


'o  laa 


1      ..      f..  .r      ^2 


Yg   .  5.  u^(H   -   Y^)      u^^^   +   G(a,T)    +  T^   Y^^^    -  TTg 
U2   =   -  I   (H   -   P)2   u^^^   +  F(a,T) 


V2   -    (H   -   P)    u^^ 


TTr,  =0 

2a 


"IT  +  %f^a   -  ^(«   -   P)^    "laaa^    +  "l   "la   + 


Substituting  for  TTp   and   p   =  Y     we   get. 


1  P 

u.  ^   +  u^   F^    -  TT-  u    (H   -   Y    )      U-,  +  U-,    U-,       +  G      + 

It  o     a        2      o^  o'         laaa  1      la  a 

u  p 

+  -^(H-Y)u,  =0 

2    ^  o  laaa 


i.e.  u,  ^+u,    u,  _+G      +uF     =0 

It  1      la  a  o  a 


^3P  =   -   "2a 


=   _    [F^   _   1    (H   -   p)2   u^^^J 
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v^  =  -  PF^  -  ^  (H  - 


')^  ^laaa  +  ^°^^^ 


0  .  -  HF^  -  1  (H  -  H)5  u^^^^  +  const 


a 


^ 


-3  =  ("  -  ^^  ^a  -  F  («  -  P)^  "laaa 


Eliminating  v^  from  this  equation  and  the  equation  (35) 


u  Y„   -  (H  -  Y  )  F^  =  -  u.  Y-,   -  Y.   - 
o   2a    ^      o   a      1   la    ix 

-  r  (H  -  Y  )^  u,     -  Y,  u, 
6  ^      o    laaa    1   la 


u   F   +  G   =  -  u,   -  U-,  Ut 
o  a    a      It    1   la 


Y 


2a 


a   2   o  ^     o    laaa    o  laaa  ~   Sa 


The  determinant 


u     Y   -  H 
o     o 


0 


u. 


o 


0 


-1 


=  0 


Hence  these  three  equations  considered  as  equations  in 
F  ,  G  ,  Y„   are  not  independent.   Adding  the  last  two 

\mAj  vX  ^_  \^ 

equations  we  have: 


Y„   +  u  F 
2a    o   a 


1  2 

-  U-,    -  U,   U-,    +  TT-  U   Ut       (H  -  Y  ) 

It    1   la    2   o   laaa  ^      o 


-  TT    +  T   Y 

Sa    o   laaa 
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Using  the  first  equations  we  have; 


u  Y„   +  u  F   = 
o  2a    o  a 


12  2 

U   U,    -  U   U,  Ut    +  tt  u   u,     (H  -  Y  ) 

o   It    o   1   la   2   o   laaa^     o 


o   Sa    o   o   laaa 


^1  ^la  -  ^IT  -  F  (^o  -  «)^  " 


-  Y   u 
laaa    1   la 


1   4 


3   v2    25 


(57)   i.e.   2Y,   =  XT  [Y.  ^  u^(i-  u   -  T  )  +  tt—   Y^  -  ^^  Y-,  +  u  tt^  ] 
^^'  '  It   da  '•  laa   o^^o    o    2u   1   u   1.   oS-' 

o       o 


Which  is  the  required  differential  equation,  determining 
the  propogation  of  disturbance.   If  we  take  the  surface 
pressure  p  and  therefore,  tTq,  constant ,  the  resulting 

O  o 

equation  is  similar  to  Korteweg-  DeVrles  equation.   The 
constant  5  can  be  explicitly  evaluated  from  the  conditions 
at  infinity.   Thus  if  the  disturbance  is  supposed  to  die  at 
a  =  -  oo ,  then  from  equations  (29),  (1?)  5  ^  *^o  '^1  *   ^^ 
now  study  time-independent  motion  with  two  special  pressure 
distributions : 


(1)      TTg  = 


(11)     TTg  = 


0      for  -  oo  <  a  <  0 
constant   for  0  <_  a  <  oo 
constant   for  |a|  <  a 
0         for  lal  >  a. 


Also  now  Y-,   =  0. 
It 


Integrating  (37)  and  using  the  fact  that  Y-j^,  ^^aa'    ^S  — ^  *^ 
as  a  — >  -  oo 
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(38) 


Y^    K  +  -^   Y? 
laa     ^.2   1 


2c. 


2c 


^1  =  -  ^S 


o 


where 


1  4 

3  o    o 


(1)    The  Integration  of  (38)  depends  upon  the  sign  of  K. 
(a)  let  K  >  0  ,  i.e.   j  c^  >  T^ 


I  (Y   -  H)^  >  T 
3^0    '      o 


1.  e , 


(n  -  d)2  >  JT. 

^  'o  gp 


i.e.  depth  of  the  water  in  equilibrium  should  be  greater 

than  s/  ~^— .   In  order  to  integrate  (38)  we  first  make  the 
gp 


substitution 
Y 


^  -  rl^i  ^o  +  T-^'^'i'^^s  '■ 


ox  2    ^ 

where  c.  >  ^  TTg.   If  c-j^  <  ^  TTg,  there  is  no 

equilibrium  position  about  which  disturbance  can  take 

place • 

The  equation  becomes 


5„  p-2  _  ^  ^c\    -   Gir^      =   0 


KG   +  -^  G 

aa    „  2      c 
2c        o 
o 


The  first  integral  is  given  by: 


(39: 


K  G^  +  ^  G^  -  ^  Ac?  -  61T, 
a    2      c   '   1     ^ 
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where  A  Is  the  constant  of  Integration. 

Let  b  =  /4c^  -  67rg  and  g  =  G  -  r,  r  a  constant 
The  equation  (59)  becomes 

K  g^  +  -4-  (g^  +  5g^r  +  ^gr*^  +  r^) 


=a   ,2 


i-  b{g^  +  2rg  +  r^)  -  A  -  0 


^o   - 

3     2 

r    br 
We  choose  r  such  that  —^  -   -  A  =  0 

c     "^o 
o 

"5        2      2 

i.e.        r-^  -  be  r   -  Ac   =  0 

o       o 


c^  \c^  0/     vc       0/ 

o  o  o 


Let  ^  =  -  3r  +  bc^      ^  =  3r^  -  2rbc 


o     -^  o   ^  o 


2    1   ,    3  ^  2  A     [Is 


o 


c^  /3K  dg 


=  a  +  B 


/-  -3  ...2 


3g   +  Ag   -  gu 
B  Is  a  constant  of  integration. 
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a  +  B  =  ± 


c  ,/3K  dg 


/--^^-- 


3g^  +  Ag   -  gu 


For  g  =  -p  , 
T 


a  +  B  =  ±  2 


Let  u  <  0 


c  /3K  dT 


{-  3  +  AT^  -  ur  ^ 


^7T72 


a  +  B  =  ±  2Cq  v/3k 


dT 


t   2c^  /3K 


1      -w( 

— 2 ^ TTir  en    - 


+A 


i^^!:i5i!i''j, 


T 


mod 


(^ 


+  A  +  /a^-12u. 
2(a2-12h)^/^ 


(40)    cn2  f  {A!zl2yJ^  (^  ^  ^^^   ^  f^  '  /^"'^^  )(G  -  r)'  . 

^2c^v/5k        3    V    ^^^    / 

Besides  B,  there  Is  another  constant.  A,  involved  Implicitly 
In  this  equation.   We  can,  however.  Impose  conditions 
say,     G(cd)  =  G'(oo)  =  0  so  that,  A  =  0.   In  that  case 
from  (39)  we  have: 
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(41) 


which  integrates  into 

^1  =  (f^o  ^1  -  f  A?  -  6Tr3  ) 

2   -,  B 

+  be   sech  i-\/bc   (— 2_  + -^l) 

^  °   c  /IT   ^ 

o' 

where  B-j^  is  again  constant  of  integration  giving 

us  one  parameter  family  of  solitary  waves.   For  some 

other  set  conditions  at  a  =  oo  for  which  !\  ^  0 , 

the  solution  is  represented  by  (40)  and  represents  the 

cnoidal  waves.   The  equation  (40)  in  terms  of  Y-,  can 

be  written  as 


^1  -ihl   'o   -f  lA^f^hr  +  ^ 

A  -  v/a^ 


12n 


^^2r(A^-12^)^/^  (g^B)} 

^  4c^/3k       J 


On  the  left  hand  side  it     =   0.   And  the  resulting 
solution  is  given  by: 

(42)       Y^  =  2c^  c-L  sech^  i   /^c^  c^  (-^—  +  B^ ) 

c/K 

where  B2  is  the  constant  of  integration  fixing  the 
position  of  a  particular  wave. 
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Let  n    >  0 


,  _  +  ^  /        dT 


a  +  B'  =  X  2 


/-  -^     ■  ^-2 


HT^  +  AT^  -  3 


For  A  <_  0,  the  above  Integral  Is  imaginary.   We,  therefore, 
assume  that  A  >  0. 


a  +  B'  =  ±  2  /   — 


/n(a2  -  t2){t2  -  p2 


Where   a?   =   ^    +  /f    "    ^^^x  p2   ^   A    -   A^    -    12^ 

1  2\i  "^1  2[i 

2 
and  A      >  12ij, .      These   restrictions    on  A   lead  to   some   restrictions 

on  A,    I.e.    on   the    Initial    conditions. 


.    =   ±  _£_  sn^   -^ 


t2 


a   +  B'    =   ±  —^ —  Sn 


°iv^        /SfT^ 


G    -   r   = 


af'   cn^  -%—   (a   +  B'  )    +  pf   Sn""  ^^  (a   +  B'  ) 


Yl   =  r   -  [-  |c^c^   +^^Ac2    -   67r3   } 


(43) 


p    p  a  [/ji  a  /nr 

^1   ""^     ~^~   («   +  B')    +  Pi   S^^   -A^   (a   +  B') 


which  represents   periodic   waves   of  period  4K  where 
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K  =  '  '  ^T 


/(I  -  t2)(1  -  "^  ^   ^  T^ 


Finally,  ij.  =  0  implies  A  =  0,  and  the  solution  in  this  case 
has  been  found  above  as  the  solitary  wave. 

11  =  5(3r^  -  2rbc^)  =  ^r-C^r  -  2bc^) 

A  =  3(-  3r  +  bc^) 

2bc 
For  M,  <  0,  0  <  r  <_ 

2bc 
But  if  r  =  — ^ —  ,  A  <  0  and  the  solution  again  is  imaginary. 

2bCo 
Hence  0  <  r  <  — ^ —  .   The  equation  giving  r  is 


(-JL.)5  _  (^L_)2  _  _A_  _ 

^bc^^     ^bc  '  ^3        ~ 

o  o  b^c 

o 

r    2 

For  ^j——   <  =-  ,   A  <_  0.   Since  the  roots  of  a  polynomial 

o    '^ 

continuously  depend  upon  its  coefficients,  A  can  be  chosen, 
i.e.  an  appropriate  set  of  initial  conditions  can  be  prescribed, 

such  that  at  least  one  root  of  the  above  cubic  will  satisfy  (44). 

2bCQ 
For  M.  >  0,  we  have  either  r  >  — = —  or  r  <  0.   But  in  this  case 

bCQ 
A  >  0,  i.e.  r  <  —^.      The  two  together  imply  that  r  <  0.   In 

2  2 

this  case  we  also  have  A   >  12ij.,  A   -  12|_l  =  9(bc  -r)(3r+bc  )  >  0 

bco 
if  '^r   +   be   >  0,  i.e.  0  >•  r  >  -  -^  .   This  again  requires 

that  A  <  0. 

The  solutions  in  the  two  cases  are  given  by  the  equations 

(4l)  and  (43)  respectively.   Each  one  of  them  involve  two 
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arbitrary  constants.   Being  periodic  solutions  these  are 

bounded  over  0  <_  a  <  oo  .   For  tTo  >  0  each  has  an  elevation 

<  2c  Ct  at  a  =  0.   This  implies  that  the  above  constants  can 
o  1 

be  evaluated  if  we  Impose  two  initial  conditions  at  a  =  0 
which  are  such  that  y-,  (0)  <  2c  c,  .   In  particular  if  we 
choose  them  to  satisfy  the  equation  (42)  we  shall  have  a 
solution  valid  on  the  whole  of  the  real  line  and  smooth  at 
least  up  to  the  first  derivative  at  all  points  of  it. 
Over  X  <_  0  we  shall  have  a  solitary  wave  and  it  will  be 
connected  on  x  _>  0  with  another  solitary  wave  or  periodic 
wave  depending  on  whether  A  =  0  or  A  ^^  0. 

(''constant   for  |a|  <_  a 
(ii)    K  >  0  TT  =   ] 

^     V.   0      for  |al  >  a 

We  again  refer  back  to  the  equation  (58)  with  ir^   as  defined 
above.   The  solution  on  the  segment  is  again  by  (41)  and 
(43)  which  are  two  parameter  families  of  periodic  waves.   On 
the  half  real  line  -  00  <  a  <  a, ,  the  solution  is  the  one 
parameter  family  of  solitary  waves 


y.  =  2c^c.  sech^  i  v^Sc  c   (-^  +  B  ) 

o'' 


As  discussed  above  these  can  be  joined  smoothly  up  to  the 
first  derivative.   We  call  the  above  solutions  W.,  ,  Wp, 
respectively.   On  the  segment  a  <  a  <  00 ,  we  have,  as  before. 
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two  parameter  family  of  waves  which  can  be,  likewise, 
connected  with  W  at  a  =  a. 

CASE  II:         K  <  0  =>  (n   -  d)^  <  -^ 

I.e.  depth  of  the  water  in  equilibrium  <  /  -^^  ,   Apparently 
this  case  does  not  arise  when  T  ==  0.   We  put  K  =  -  ^, 
i  >  0  In  the  equation  (38) 


-z  p    2c -, 

laa   2c2   1    %   1    S 


(1)  TTg  =  0   for   -00  <  a  <  0 


so  that 


2^0         ° 
We  write  Y^  =  -  Z-j^.   The  above  equation  becomes 


2c^       o 

This  equation  represents  solitary  wave  only  If  c,  <  0. 
Then  It  Integrates  Into, 


Z  =:  -  2c.  c 


2  1  ./-^ ,  a 


1   o 


ech-±/r2F7T;  (-^+d^) 


(46) 

Y-,    -   2Ct    c^    sech^  J  i/-2c,    c       (- 


1   -   2c^   c^    sech^   2-  1/-2C1    c^    (-^  4-   d^) 
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Since  c-,  <  0,  this  Is  a  negative  wave.   d..  Is  constant  of 
integration  which  would  represent  a  parameter. 

On  0  <  a  <  00  ,  ir^  j^  0.      The  solution  as  before  is 

-  ^1  =  -  I  ''l  ''o  -  T  /^4  -  ^"^S  +  ^  +  — ^ 


/^- 


12u 


^2   (^"-l^M-)'^"^  (a  +  B) 


4c^  /3l 


c     ,- 
(47)  ^i   =  ^  ^i   %   +   -^v/^^1  -  ^Trg  +  r  +  ^ 


A  -  /a^-12^1. 


^2AAf-12tU^( 


where  A,  |j,,  r  have  the  same  meaning  as  before.  Therefore 
the  result  in  this  case  also  is  that  the  two  waves  can  be 
joined  smoothly  at  a  =  0.   Similarly  in  (11)  when 

constant   for   |a|  <  a 
0       for   I a|  >  a 

We  get  the  three  solutions  that  can  be  connected  smoothly 
up  to  the  first  derivative. 

Finally,  the  complete  solution  consists  of  two  or  more 
solitary  waves  or  the  solitary  waves  and  the  periodic 
waves  joined  smoothly  up  to  the  first  derivative  at  the 
points  of  discontinuities  of  the  pressure.   The  levels  at 
the  two  infinities  need  not  be  the  same.   The  solitary  waves 
can  be  positive  or  negative  depending  on  whether  the  depth 
at  a  =  -oo  is  greater  than  or  less  than  \/  ^  . 

V   gp 
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Chapter  III 
Waves  In  Shallow  Waters 

Formulations  of  the  problem 

We  consider  water  flowing  In  an  Infinite  channel 
with  finite  depth;  the  bottom  of  the  channel  Is  considered 
horizontal.   The  assumptions  on  the  flow  and  the  water 
are  the  same  as  Imposed  analogously  In  Chapters  I  and  JI. 
We  take  the  x-axls  In  the  bottom  plane  and  y-axls 
vertically  upward.   These  axes  are  considered  fixed  In 
space.   Let  the  free  surface  be  given  by  y  -  t^(x,  t). 
The  Euler   equations  of  motion  are: 

(1)  u,  +uu+vu=-  —  p 

^    '  t     X     y     p  ^x 

(2)  V,  +uv+vv= p-g 

^  t     X     y     p  y   ^ 

The  equation  of  continuity  is 


(3)  ^x  +  V  ^  ° 


Kinematlcal  condition  at  the  free  surface  is 


(4)  ^^  +  UT]^  -  V  =  0 

V  =  0  at  bottom. 

(5)  ^  -  Po=    ,,l   g'^yg       y  =  ^^ 
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where  p^  is  the  atmospheric  pressure  at  a  point  just 
outside  the  liquid  and  ir  Is  the  pressure  at  the  correspond- 
ing point  just  Inside  the  liquid. 

Integrating  the  equations  of  continuity: 


V 


u  dy 


0 


Substituting  for  v  in  (k) : 


0 


+  UT] 


X 


u^  dy 


0 


Also  we  have: 


-^   /  u(x,y,t)  dy  = 

b 


"x  ^^   ^  Tlx^(x,y,t) 


0 


y=Ti 


The  two  together  =- 


^t  =  ^  /  ^  ^y 


0 


For  the  shallow  wave  theory  we  have  the  hydrostatic  law: 


P  -  Po  =  - 


TT] 

2^%^  +  gp(Ti-y) 


It  is  obvious  from  this  that  p   does  not  depend  upon  y. 
Hence  If  u  is  originally  Independent  of  y,  it  will  be  so 
always.   We  suppose  it  likewise  so  that  u  =0. 
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The  equations  then  reduce  to 


i.e. 


1  B   /     ^xx 

^t  +  ^^x  +  s^x  =    ?  ^  ^TT^Tp/^' 


(6)  U^  +  UU^  +  gTl^  =    -  ^  (^ 27^75) 

'x' 


(7)  -\=^^-^^ 


If  we  suppose  that  water  has  zero  surface  tension,  then 
It  Is  well  known  that  the  equations  (6),    (?)  with  t  =  0 
has  uniform  flows  as  Its  only  smooth  solution.   The  problem 
dealt  with  here  Is  the  investigation  of  some  other  smooth 
flows  when  t  7^  0 .   In  particular  we  look  for  progressive 
wave  solutions. 
Let 

•n  =  Ti(x  -  at)  =  Ti(0 

u  =  u(x  -  at)  =  u(C) 

The  above   equations    then  become. 


(1^^) 


ari '    =  -^    {uT]) 
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Integrating,    we  get 


(8)  ^    (u-  a)^   +  gTi   =   +  ^     p%5  +  A 


T)  ^^ 
P  (1+Tl|)" 


27372 


(9)  at)   =  UT]    +  B 

We  assume  that  at  C  =  0^  "H  =  h,  tj  '  =  b,  !]"=  0,  u  =  u^, 
b  is  some  positive  constant.   Then 

A  =  2  (Uq  -  a)   +gh 
B  =  (a  -  UQ)h  . 
Eliminating  (u^-a)  from  (8),  (9)  we  get 

Let  T]  '  =  V;  the  result  is  a  first  order  ordinary 
differential  equation: 

p   dT]       ^^2     ^^'     2 
(10) 

T   dV  _  (1+V^)^/^  ,         2   1   ,     .2 
p   d7[  -  ^^~2     (h-Ti)[-gTi  +  2  hlUQ-a) 

+  -^  T](uQ-a)  ]  . 

This  equation  governs  the  disturbance  about  the  following 
equilibrium  positions  which  are  the  singular  solutions -.V  =  0, 
T)  =  h  and 


7^ 


-  I  (^o-a)^  +  {  ^  (uQ-a)^+  2gh(uQ-a)2]l/2 

-  2g 


4g 


^4-'  n. 


^ 


where  ti  =  •q   corresponds  to  the  positive  sign  and 

Ti  =  T]_  corresponds  to  negative  sign. 
The  equilibrium  positions,  therefore,  are 


V  =  0,      T]  =  h,  Tl_^,  T]_ 

To  start  with  we  consider  the  disturbance  about  t]  =  h. 
Set  A  =  ri-h. 

T   dV      (1+V^)^^^A   r   /.^v.^2^  1  v,^     ^2 
-  TTT  =  -  -^ —^ [-g(A+h)  +  p  h(u^-a) 

+  I  (uQ-a)^(A+h)} 
(11) 

V(A2+2Ah+h^) 


[ -gA^+A^{ -2gh+  ^(u^-a) 2] +A{ h(u^-a) ^-gh^} ] 


V(A^+2Ah+h^) 


from  which  it  is  obvious  that  in  order  to  consider  the 
behaviour  of  the  solution  about  A  =  0  we  have  to  study 
the  equation: 
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^  ^  =  -  Ah{(uQ-a)^-gh}/Vh^. 
We  consider  the  two  cases: 

o 

(a)  (uq  -a)  ^   gh 

(b)  (uq-  a)   =  gh  . 

2 
Case  a.   For  (up^-a)   >  gh,   t]  =  h   Is  a  center. 


^0 

^0- 


2 
For  (up^-a)   <  gh,   t]  =  h  is  a  saddle  point 


2 
For  (upj-a)   <  gh. 


(up^-a)   +  [  (up,-a)   +  8  gh(up,-a)  }  ' 


Tl 


+  -  irg 

V,  ^  (    2,2^p  2,  2  ^  1/2 
^  gh  +  (g  h  +8g  h  )      ^  ^ 

^fg 

Therefore  ri  ,  <  h.   Also  ri   <  0.   Hence  r]   <  t]   <  h, 
+  -  —    + 

For  (uQ-a)^  >  gh. 


Tl 


(uQ-a)2  +  {(uQ-a)^  +  8gh}l/2 

>  P^\+  ^^^  ^  h. 
4g 


In  this  case  t]_  <  h  <  t]  . 

We  now  consider  the  singularity  at  t]  =  tj  . 


^f  ^-Il^^(.-,)(n-..)(.-V 


Then 


78 


g{(l+V^)^/2-l}{A(A+A_^-h)(A+n^-Ti_)} 

gA(A+Tl      -h)       (A+Tl      -T]      ) 

4-  21 Z I ^ 

Whence  it  follows  that  we  have  to  consider  the  linear 
equation: 


o 
Implies  T)  =  Ti   Is  a  saddle  point  for  (u^^-a)   >  gh 

Similarly  for  ri  =  ri_  we  consider  the  linearized 
differential  equation: 


p 
and     T  =  "H   Is  a  center  for   (up^-a)   <  gh. 


gA  {t]_    -   h)  (ti_  -  T)^) 


T  dV  ^  

P   dA   -  ^  ^2 

I.e.  Tj  =  ri_  Is  always  a  saddle  point. 
We  know  from  Polncare's  theory,  that  a  center  In  the 
linearized  differential  equation  goes  Into  a  center  or 
spiral  In  the  nonlinear  case.   It  Is,  therefore, 
necessary  to  further  analyze  the  differential  equation 
(11)  to  determine  the  nature  of  Its  singularities  which 
are  analogous  to  centers  In  the  corresponding  linearized 
differential  equation  (12) . 
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For  this,  we  integrate  (10), 


'^  V       dV  _       1  B^    „ 

P   (1+7^)5/2   d,  -  S^  +  2  ^2  -  ^ 


P    /— 7:2        2     2  "1 


2       2 
=  -^3 i ATI  +  const 


1+  V 


But   C   =  0,    ri   =   h,    T]'    =  b,  u  =  Uq.      Therefore 

-  ^     1 =   gh!  1   B^   _   Ah   +   const 

P        / o                2  2    h 

yi  +  ^ 


P      / p  / 2  ^         Ti        n 

/  1+b^         /l   +V^  2 

-,                                  h(up^-a) 
=  -|(Ti-h)[g(Ti+h)+  ^ 

p 
-    (uQ-a)      -   2gh] 


I.e. 

We  write  gL   =    (uQ-a)      -  gh. 


/T+b^  /  1+V 


~^]    =  ^  -^^i^    (gri-gL-gh) 
2  ^  I 


Let  a    =    pg/2T; 


-^  =  ^  -  ^iH^   [  (n-h)-L}    ; 

2  T 


/l+V^       /  1+b^ 
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1  +  v^    =   n 


[  — zL=  -  a(Ti-h)2(Ti-h-L)]2 

In  order  to  see  whether  the  singularities  r]  =  h,  r) 
are  centers  or  spirals.  In  the  nonlinear  case,  we 
have  only  to  take  the  Intersection  of  this  curve  with 
a  line  passing  through  t^  =  h,  or  t]  =  ti_|_  as  the  case 
may  be;  I.e.,  with  V  =  m(ii-h),  V  =  m(T]-T]_^)  where  m  Is 
a  parameter  which  rotates  the  line  through  2v. 
Thus  we  have: 


i+m^(Ti-h)^    = . n 


[    ''    -  a(Ti-h)^(Ti-h-L)]^ 

which  for  any  fixed  m.  Is  an  eighth  degree  In  t]   and 
so  can't  have  more  than  eight  roots.   The  spiral,  on 
the  other  hand,  must  be  Intersected  by  any  radial  line 
In  Infinitely  many  points.   Hence  In  both  the  cases  the 
singularities  at  t]  =  h,  t]  =  i^^  must  be  centers.   See 
Figures  1  and  2  below.   The  saddle  point  shows  that  the 
disturbance  Initiated  at  t]  =  h  reaches  this  level  again 
at  ^— >-oo  or  ^— *--oo,  depending  upon  the  sign  of  b. 
Thus  In  our  case.  It  tends  to  form  a  crest  on  the  right 
and  a  hollow  on  the  left  of  the  x-axls .   In  each  case 
there  Is  a  single  hump  or  hollow. 
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Case  b.    (up  -  a)   =  gh. 

The  differential  equation  In  this  case  becomes: 

T   dV         (1+V^)^/^    f^      \/         2^  1      ,  2    ,    1      ,     ^ 
■p   dTi'  "  "^ 5 (h-Ti)(-gTi    +  2  gh     +  -^  ghTi) 

^   g(l+V    ) (h-Ti)(-2ri^    +  h^    +      h-q) 

2Vt]'^ 

^  _£U+V) (h-Ti)(h-Ti)(2Ti+h) 

2Vti 

(11)  f    =    f  ■    iiJ:!^  (n-h)2(2,+h). 


The  singular  points  are  V  =  0,  i]  =  h  and  V  =  0,  t)  =  -h/2 , 
The  former  Is  a  double  singular  point.  Integrating  (11), 
we  get  ■ 

l+V^  =  IL 


[     ^     -  a(Ti-h)5]2 
2 


/  1+b' 


where  a  =  pg/2T. 

In  order  to  discuss  these  curves  In  the  vicinity  of  r)  =  h, 
we  have  the  substitution  rj-h  =  A. 
Then  we  have: 
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1  +  v^  =  ^^+^)^ 


[  — — --  -  aA^] 


J^ 


I.e. 

1  +  V^  .       1+^^ 


^■^     A+h   '^  ■■ 

where  A  Is  considered  small. 

This  then  in  the  vicinity  of  A  =  h,  yields; 


1   +  V^  ==    (l+b2)    [1   +  2a^^3] 


h 


V2   =   b2    +    2^    (l+b2)5/2   ^3 

h     ^  ' 

2  _  2a    m^h2,3/2  (- hb^ ^_  ^  ,3, 


V"   = 


^  ^2a(l+b2)^/2 


d^./f   U+b2)5A    [,3,,3p/2 


5  .  hb^ ^^^     R     _   /^  n..K20/4 


where     ^i''  =  5-^T^   •  Let     B,    =     ^   (l+b    ) 

Then 

Integrating,  we  get 


B-,^  +  const.  =  / 


A 

dA 


Ju^TT^ 
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We  have   to  evaluate   the   Integral: 


r 


d?\ 


dA 


A^+U^  /(A4u)      /  A^-Au  +  u^ 


-2 


dT 


yF^-  -  .2 


'3uT  +1-3uT' 


A+M-   =  — 7j 
T 


dT 


[ (T^-  y^rT+r) (T^+  y3rT+r ) ] ^ ^^   ' 


A^ 


'3H' 


j^  r 


dT. 


2+  /3         -^         2-/3 


1/2 


dT. 


,   ^   T+  A^ 
'^      T-  /F 


v'^     ^        {(T^+(2-   ,^)2)(t2+(2+/3)^)1^^^ 


k  r 

/ 


dS 


./?(2+  /3) 


yw  /i--i^ 


(2+/3)" 

( 

for    (2+/5)^+T^   =    ^^'^{^^' 


4  _-] 

B-,4    +   const.    =  Sn      S 

^  v^    (2+/3) 


Sn 


-1  (2+y3)(T-  y^) 


?    (2+,/5)  {(2+/3)^(T-y?)2+(T+yr)^}^/^ 
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where  d  is  an  integration  constant. 

(TH-  /?)^      ^  cn^t  -^    (2+/3)(p^g  ^-  d)) 


=> 


r(-n-h-Hj.)  =  [ — _ -jS 

where  r  =  — ^  ,   p.^  =  ^^  ^  ,  .„  ^  ct  =  -gf 

3,=  /f  (l+b2)5A 


=>  that  we  have  periodic  waves  of  period 

C2+  /3)(2a)^/5 

1 

where  ^  =    f  ^^         ,   k^  =  ^  /3 

0   y(l-s'^)(l-kV)  (2+/5)2 

The  other  singularity  is  at  V  =  0,  ri  =  -  b/2. 
This  is  rejected  in  view  of  the  fact  that  this  is 
below  the  bottom  of  the  channel. 
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Thus  we  have  the  following  result: 
For  T  =  0,  we  have  uniform  flows  as  the  only  smooth 
flows.   When  x  ^  0,    we  divide  the  solution  In  three 
parts.   For  (u^  -  a)  <  gh,  the  disturbance  starting  at 
a  point  Is  carried  over  on  both  the  sides  of  the  source 
of  disturbance.   This  disturbance  propagates  In  a 
single  crest  on  one  side  and  in  a  single  trough  on  the 
other  side  right  up  to  +oo  or  -co  as  the  case  may  be. 
In  case  (u^^  -  a)   >_  gh,  the  disturbance  is  propagated 
as  a  periodic  wave. 
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13  ABSTRACT 


The  purpose  of  the  present  work  Is  to  study  surface  waves  in 
rectangular  channels  of  finite  depth.   In  all  cases  we  consider  two 
dimensional  flow  of  an  Incompressible  and  nonvlscous  fluid.   We  also 
suppose  that  there  is  surface  tension  present  in  it. 

In  the  first  chapter  we  discuss  the  feasibility  of  a  steady 
flow  consistent  with  the  linearized  problem,  and  study  the  corre- 
sponding radiation  conditions.   Such  a  flow  is  found  to  be  possible 
In  both  the  supercritical  and  the  subcritlcal  cases  but  not  in  the 
critical  case.   In  the  supercritical  case  the  disturbance  dies  both 
upstream  and  downstream;  in  the  subcritlcal  case,  in  general.  It 
does  not  vanish  downstream. 

The  second  chapter  deals  with  a  study  of  solitary  and  periodic 
waves  which  result  when  a  stream  of  uniform  flow  Is  disturbed  by  an 
arbitrary  variable  pressure  on  its  free  surface.   We  obtain  a 
general  differential  equation  of  the  motion.   Explicit  wave  solutions 
are  obtained  by  considering  two  special  pressure  distributions. 

In  the  final  chapter  we  study  the  progressive  wave  solutions  of 
the  zero  order  shallow  water  equations.   In  the  subcritlcal  case,  we 
obtain  a  wave  with  a  single  hump  and  a  single  hollow  on  the  two 
sides  of  the  point  of  disturbance.   In  the  critical  and  the  super- 
critical cases,  we  obtain  periodic  waves. 
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